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Preface 


This  volume  contains  Appendices  A  through  D  that  were  referenced  in  Volume  1; 
Technical  Summary  of  the  ClAO  II  workshop  results.  Appendix  A  provides  overviews  of 
the  ClAO  I  and  ClAO  II  workshops,  participants  and  summaries  of  papers  presented  at 
the  ClAO  I  workshop.  Appendix  B  provides  an  extensive  bibliography  of  millimeter  wave 
data  reports  for  laboratory,  field,  and  operational  measurements.  Appendix  C  compares 
results  from  the  model  comparison  exercise  conducted  for  the  ClAO  II  workshop. 
Appendix  D  provides  a  review  of  existing  measurement  techniques  and  facilities  where 
these  techniques  are  employed. 
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1.  Purpose  and  Scope 


Appendix  A  provides  overviews  of  the  CIAO  I  and  CIAO  II  workshops.  A  summary  abstract  of 
the  work  discussed  and  problems  posed  by  the  participants  is  provided. 

First,  CIAO  I  and  CIAO  II  participants  and  their  mailing  addresses  are  given.  Next,  presentation 
overviews  from  CIAO  I  that  summarize  the  thrust  of  the  research  discussed  during  this  workshop 
follow.  CIAO  II  did  not  use  a  paper  presentation  format  and  results  from  the  discussions  held 
during  that  workshop  are  described  in  Volume  1  of  this  report.  The  CIAO  I  presentation 
overviews  are  organized  into  three  areas:  1)  model  development,  2)  model  application,  and  3) 
experimental  evaluation.  Presentations  on  model  development  discuss  the  development  of 
models  of  scattering  and  absorption  from  first  principle  physical  assumptions.  Papers  on  model 
applications  discuss  the  application  of  previously  developed  models  to  the  analyses  of 
experimental  data  or  particular  problems.  Papers  on  experimental  evaluation  discuss 
experimental  methods  and  analyses. 
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1.3  CIAO  I  Presentation  Overviews 

The  CIAO  I  presentation  overviews  are  organized  into  three  areas:  1)  model  development,  2) 
model  application,  and  3)  experimental  evaluation.  Papers  on  model  development  discuss  the 
development  of  models  for  scattering  and  absorption  from  first  principle  physical  assumptions. 
Papers  on  model  applications  discuss  the  application  of  previously  developed  models  to  the 
analyses  of  experimental  data  or  particular  problems.  Papers  on  experimental  evaluation 
discuss  experimental  methods  and  analyses. 

1.3.1  Model  Presentation  Overviews 

Overviews  presented  in  this  section  represent  those  describing  model  development.  In 
overviews  representing  the  work  of  more  than  one  author,  the  workshop  presenter  is  underlined. 
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ADVANCED  COMPUTATIONAL  TECHNIQUES 

by 

Dr.  John  Fratamico 

Science  Applications  International  Corporation 
San  Diego,  CA 

Dr.  Fratamico’s  experience  is  in  the  design  and  development  of  Stealth  weapons 
platforms.  His  talk  was  focused  on  software  and  hardware  optimization  techniques  developed 
in  stealth  technology  that  may  have  application  to  atmospheric  physics  problems.  The  specific 
scattering  technique  referenced  is  the  three  dimensional  method  of  moments  (3D-MoM).  This 
is  a  technique  that  spans  the  full  gamut  of  material  properties,  material  complexities  and  shapes. 
Advanced  computational  techniques  addressed  effects  that  can  broaden  the  particle  size  and 
range  to  encompass  as  broad  a  range  of  material  parameters  as  possible. 

The  first  topic  considered  was  the  antenna  type  problem  where  the  wavelength  is  large 
compared  to  the  particle  size.  For  this  case,  the  impedance  boundary  condition  formulation  is 
appropriate  as  one  may  vary  not  only  material  properties  but  also  the  impedance  loading. 
Computationally,  the  formulation  should  avoid  _-n’B  terms  when  the  particle  size  is  much  smaller 
that  a  wavelength:  avoid  B-n’B  terms;  and  force  the  impedance  terms  near  the  diagonal  of  the 
matrix. 


A  typical  out-of-core  impedance  matrix  might  consist  of  100,000  x  100,000  fully  complex 
elements.  This  requires  approximately  100  gigabytes  of  storage.  The  MoM  solution  requires 
on  the  order  of  floating  point  operations  and  I/O  operations.  Efficient  algorithms  for 
interleaving  these  operations  on  both  parallel  and  vector  architecture’s  were  discussed. 

Another  important  point  addressed  was  the  limits  of  excursion  from  an  exactly  solved 
problem  that  may  be  pursued  via  perturbation  techniques.  If  the  problem  can  be  reconfigured 
such  that  the  known  problem  is  pushed  into  the  lower  right  corner  of  the  MoM  matrix,  then  one 
can  obtain  the  scattered  fields  for  the  perturbed  problem  without  solving  the  intermediate  step 
for  the  surface  currents.  This  reduces  computer  run  time  by  a  factor  of  approximately  1/1,000. 
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THE  DISCRETE  DIPOLE  APPROXIMATION  FOR  SCATTERING  CALCULATIONS 

by 


Bruce  T.  Draine 

Princeton  University  Observatory 
and 

Piotr  J.  Flatau 

Scripps  Institution  of  Oceanography 
March  1994 

The  discrete  dipole  approximation  (DDA)  is  an  approximation  of  a  continuum  target  by 
a  finite  array  of  poiarizable  points.  The  points  acquire  dipole  moments  in  response  to  the  local 
electric  field.  An  array  of  polarizable  points  can  accurately  approximate  the  response  of  a 
continuum  target  on  length  scales  large  compared  to  the  dipole  separation.  For  a  finite  array 
of  point  dipoles,  the  scattering  problem  can  be  solved  essentially  exactly.  The  only 
"approximation"  in  the  DDA  is  the  replacement  of  the  continuum  target  by  an  array  of  N  point 
dipoles.  The  replacement  requires  specification  of  both  the  vector  locations  of  the  dipoles  and 
the  dipole  polarizabilities.  There  is  a  close  correspondence  between  the  DDA  and  discretizations 
based  on  the  "digitized  Green’s  function"  method  or  the  "volume  integral  equation"  formulation. 
In  the  low  frequency  limit  where  the  interdipole  spacing  or  side  of  a  cubical  subvolume  is  much 
less  than  an  incident  wavelength  the  DDA  has  been  shown  to  be  formally  equivalent  to  digitized 
Green’s  function  and  volume  integral  formulations.  However,  workers  developing  these 
formulations  encountered  difficulties  in  attempting  to  compute  second  order  corrections  whereas 
DDA  permits  rigorous  derivation  of  second  order  corrections.  The  principle  limitation  of  the  DDA 
is  in  handling  the  target  boundaries  since  the  DDA  array  has  a  minimum  length  scale  equal  to 
the  interdipole  spacing.  For  targets  with  large  values  of  the  refractive  index,  the  accuracy  of  the 
DDA  suffers.  In  this  case  the  DDA  can  seriously  overestimate  the  absorption  cross-section,  even 
in  the  dc.  limit.  When  the  absolute  value  of  the  refractive  index  is  much  greater  than  one  other 
techniques  may  be  superior  to  the  DDA. 

Fast  Fourier  transform  (FFT)  methods  can  greatly  accelerate  the  computations  required 
to  solve  the  scattering  problem,  but  only  if  the  dipoles  are  located  on  a  periodic  lattice.  With  FFT 
methods  calculations  can  be  performed  on  a  workstation  for  numbers  of  dipoles  (e.g.  ~  10®) 
much  greater  than  the  largest  values  (~  1 0^)  which  could  be  handled  without  FFT  methods.  We 
have  used  FFT  methods  accepting  the  requirement  of  a  periodic  lattice  and  further  assumed 
cubic  lattices. 

There  have  been  questions  of  the  best  method  for  assigning  dipole  polarizabilities.  One 
approach  is  to  use  the  "Claussius-Mossotti"  polarizabilities  which  is  exact  for  an  infinite  cubic 
lattice  in  the  limit  that  the  interdipole  spacing  is  much  less  than  a  wavelength. 

Extensive  DDA  calculations  for  spheres  comparing  different  prescriptions  for  the  dipole 
polarizability  confirms  that  the  lattice  dispersion  relation  developed  analytically  by  Draine  and 
Goodman  for  an  infinite  lattice  of  polarizable  points  appears  to  be  the  best  when  the  product  of 
the  absolute  value  of  the  index  of  refraction,  wave  vector,  and  interdipole  spacing  is  less  than  1 . 
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By  locating  dipoles  on  a  lattice,  computational  memory  requirements  are  approximately 
linearly  dependent  on  N,  the  number  of  dipoles.  For  example,  the  program  DDSCAT  has  a  total 
memory  requirement  of  approximately  O.SSfN^NyN^/IOOO)  Mbytes,  where  x  Ny  x  is  the 
rectangular  volume  containing  ail  the  N  dipoles.  CPU  requirements  are  also  significant  for  large 
targets.  For  example,  on  a  Sun  4/50  (Sparcstation  IPX),  a  single  CCG  iteration  requires 
appro.ximately  3.0(N,NyN,/1000)  cpu  seconds.  Between  10  and  100  iterations  are  typically 
required  to  solve  for  a  single  incident  direction  and  polarization.  It  is  for  this  reason  that  T-matrix 
methods  which  exploit  efficient  procedures  for  orientational  averaging  are  competitive  for  some 
target  geometries  as  well  as  recursive  T-matrix  algorithms  currently  developed  by  Chew. 

At  this  time  there  is  no  known  way  to  predict  precisely  the  accuracy  of  DDA  calculations. 
Computations  for  the  extinction  efficiency  with  moderately  refracting  spheres  performed  using 
DDA  and  compared  with  exact  Mie  theory  calculations  show  errors  are  less  than  3%  for  size 
parameters  less  than  12  and  an  index  of  refraction  of  1.33  +0.01  i.  For  an  index  of  refraction  of 
2+i,  the  error  in  the  absorption  efficiency  is  about  4%  when  dipole  numbers  are  greater  than  1 0“. 
Computations  for  the  error  in  the  differential  scattering  cross-sections  are  within  20%  for  all 
scattering  directions  for  an  index  of  refraction  of  1.33+0. 01  i  and  within  30%  for  2+i. 

The  DDA  has  been  applied  to  calculations  of  scattering  and  absorption  by  rough  and 
porous  particles,  interstellar  graphite  particles,  aggregate  particles,  inhomogeneous  particles, 
structures  on  surfaces,  microwave  scattering  by  ice  crystals,  single  scattering  properties  of  ice 
crystals,  and  scattering  at  plane  boundaries  perturbed  by  roughness. 
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DIGITIZED  GREEN’S  FUNCTION  APPROACH 
TO  SCATTERING  FROM  ARBITRARY  STRUCTURES 

by 

Dr.  George  Geodecke 
New  Mexico  State  University 
Las  Cruces,  NM 

The  digitized  Green’s  function  approach  to  scattering  from  complex  objects  was 
discussed.  The  method  is  completely  general  and,  in  principle,  is  capable  of  calculating 
absorption  and  scattering  by  waves  of  all  kinds  from  localized  scatterers  of  arbitrary  shape,  size 
and  morphology.  Mathematically,  the  digitized  Green’s  function  approach  is  closely  related  to 
the  coupled  dipole  method  (Purcell  and  Pennypacker,  AP  J  1^,  705-714  1973). 

The  development  discussed  at  this  conference  included  the  proper  way  to  include  self¬ 
terms  and  depolarization  dyadics  into  the  formalism.  It  was  shown  that  the  solutions 
automatically  obey  the  Optical  Theorem.  Two  specific  applications  were  discussed:  1 )  structures 
with  smoothly  varying  constants;  and  2)  structures  with  a  few  sharp,  cutoff  surfaces  such  as  soot 
or  snowflakes. 

Areas  identified  as  needing  further  research  included  incorporating  a  Born  series 
approximation,  better  accounting  for  near  neighbor  interactions,  and  using  a  variable  cell  size 
over  the  parameter  space.  Limitations  of  the  technique  were  also  discussed.  These  are 
basically  related  to  computation  time  and  storage  limitations.  These  issues  were  also  addressed 
by  Dr.  John  Fratamico. 
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SIMPLIFIED  TREATMENT  OF  THE  EXTENDED  QUASISTATIC  THEORY 
FOR  ABSORPTION  BY  CONDUCTING  FIBERS 

by 

Norman  Pedersen,  Jeanne  Pedersen, 
and  Peter  Waterman 
Panametrics,  Inc. 

26  June  1989 

An  elementary  treatment  of  electromagnetic  wave  absorption  by  conducting  fibers  is 
developed  assuming  the  fibers  are  prolate  spheroids,  the  quasistatic  approximation,  that 
wavelengths  are  greater  than  20//m.  We  note  that  even  in  the  low  frequency  limit,  the  problem 
of  analytically  calculating  the  electric  field  distribution  within  a  right  circular  cylinder  has  not  been 
solved  in  closed  form.  The  quasistatic  approximation  requires  the  cylinder  length  to  be  much 
less  than  an  incident  wavelength.  However,  in  spite  of  the  assumptions  made  in  the  analysis, 
accurate  results  for  absorption  by  very  thin  fibers  are  obtained  over  very  wide  ranges  of 
meaningful  wavelengths  and  easily  obtained  particle  dimensions.  In  many  cases,  one  or  more 
of  the  limiting  assumptions  can  be  relaxed  without  adverse  results. 

For  randomly  oriented  fibers  of  high  aspect  ratio,  the  absorption  cross-section  has  the 
following  features: 

1 .  When  the  product  of  the  depoiarization  factor  and  imaginary  component  of  the 
dieiectric  constant  squared  is  much  less  than  1,  large  values  of  particle  material  electrical 
conductivity  produce  large  values  of  the  absorption  cross-section  per  unit  volume. 

2.  Conditions  for  feature  1,  can  be  achieved  by  making  particle  aspect  ratio  large. 

3.  High  values  of  absorption  cross-section  per  unit  volume  obtained  using  features  1  and 
2  are  the  result  of  a  non-resonant  condition  that  can  be  met  over  a  wide  range  of  frequencies. 

The  quasistatic  approximation  has  been  extended  to  provide  an  analysis  of  the  absorption 
cross-section  when  the  particle  length  approaches  or  exceeds  the  free  space  wavelength  of  the 
incident  radiation.  The  extended  quasistatic  approximation  assumes  the  effective  cylinder  half- 
length  is  half  a  wavelength.  This  assumption  cannot  be  supported  on  strict  theoretical  grounds 
but  has  been  found  to  provide  accurate  and  useful  results  in  many  cases.  The  result  of  our 
analysis  with  this  approximation  is  to  slightly  modify  the  depolarization  factor  used  in  quasistatic 
approximation  computations.  The  advantage  of  the  Extended  Quasistatic  result  is  that  it  allows 
one  to  use  analytic  methods  in  obtaining  closed  form  approximate  solutions  to  practical 
problems. 
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INTEGRAL  EQUATION  ANALYSES  OF  ARTIFICIAL  DIELECTRICS 

by 

Dr.  Edward  Newman 
Ohio  State  University 
Columbus,  OH 

Tl^e  electromagnetic  analysis  of  artificial  dielectrics  is  one  of  the  fundamental  building  blocks 
of  Low  Observable  (Stealth)  Technology.  Dr.  Newman  and  his  group  at  Ohio  State  University 
are  the  premier  leaders  in  the  unclassified  research  into  this  area.  The  techniques  discussed 
are  important  and  directly  applicable  to  problems  in  atmospheric  physics.  The  electromagnetic 
analysis  of  an  artificial  dielectric  and  the  electromagnetic  analysis  of  a  heterogeneous 
smoke/obscurant  cloud  are  intimately  intertwined. 

This  paper  focuses  on  two  mathematical  techniques:  1)  the  formulation  of  the  relevant 
integral  equations,  and  2)  their  solution  by  Method  of  Moments  (MoM)  techniques.  Calculated 
results  are  shown  to  compare  favorably  with  measurements  for  three  cases:  1)  a  fixed  array  of 
dipoles,  2)  a  fixed  array  of  crosses,  and  a  3)  composite  weave.  These  are  directly  applicable 
to  the  Low  Observables  problems  of  conformal  antennas,  band-pass  radomes,  and  wing  leading 
edge  treatments,  respectively.  The  extension  to  the  complex  artificial  dielectric  properties  of  a 
heterogeneous  aerosol/particulate  cloud  is  apparent. 

The  examples  discussed  in  this  paper  address  the  effects  on  the  effective  permittivity  due 
to  the  host  medium,  the  scattering  medium,  the  frequency,  the  effective  lattice  spacing  of  the 
artificial  dielectric,  the  scatterer  shape  and  the  polarization  and  angle  of  the  incident  wave.  For 
the  undriven  case,  it  is  shown  that  the  resulting  integral  equations  are  homogeneous  and  the  root 
corresponds  to  a  different  polarization  of  the  incident  wave. 

It  is  demonstrated  by  comparison  with  experimental  results  that  the  techniques  described 
accurately  account  for  complex  shapes  of  the  individual  scatterers  and  the  mutual  coupling 
between  them.  It  is  explicitly  shown  that  the  mutual  coupling  between  individual  particles  is 
adequately  accounted  for  by  increasing  the  number  of  expansion  functions  used  on  the  MoM 
solution  to  the  integral  equations. 
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WAVELET  TRANSFORMS  IN  THE  T-MATRIX  THEORY 
OF  ELECTROMAGNETIC  SCATTERING 

by 

Dr.  Patrick  L.  Nash 
University  of  Texas  at  San  Antonio 
San  Antonio,  TX 

The  construction  of  the  T-matrix  requires  the  numerical  evaluation  of  many  surface 
integrals  involving  the  tensor  products  of  the  vector  spherical  harmonics  and  the  angular 
dependent  factors  describing  the  boundaries  of  the  scattering  particles.  The  wavelet  transform 
approximates  these  square  integrable  functions  with  smoothed  versions  of  these  functions.  This 
results  in  a  matrix  formulation  that  lends  itself  to  a  multi-wavelength  resolution  analysis  of  the 
scattering  problem  under  consideration. 

This  paper  presents  results  obtained  from  using  a  generalized  two-dimensional  wavelet 
transform  and  an  approximate  matrix  inversion  technique.  Iterations  of  the  wavelet  transform 
on  an  n  X  n  matrix  are  shown  to  yield  a  sparse  and  easily  solvable  problem  even  in  the  case 
where  the  initial  diagonal  elements  are  arbitrary  (and  even  singular). 
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LIGHT  SCATTERING  FROM  AGGLOMERATES: 

COUPLED  ELECTRIC  AND  MAGNETIC  DIPOLE  METHOD 

by 

George  W.  Mulholland 

National  Institute  of  Standards  and  Technology  (NIST) 

Gaithersburg,  MD 

This  is  an  outgrowth  of  a  fire  research  program  at  NIST.  Of  particular  interest  is  the 
growth  of  soot  and  smoke  particulate.  It  is  shown  that  the  absorption  cross  section  of  a  complex, 
fractal  soot  particulate  is  approximately  given  by  the  absorption  cross  section  of  the  basic 
building-block  structure  times  the  number  of  these  building  blocks  in  the  agglomerate.  It  is 
usually  sufficient  to  use  a  sphere  as  the  basic  building  block  and  to  model  the  agglomerate  as 
a  collection  of  touching  spheres.  No  multiple  scattering  within  the  agglomerate  is  considered. 

Inclusion  of  a  magnetic  dipole  term  in  the  solution  is  shown  to  essentially  double  the 
range  of  size  parameter  for  which  the  technique  yields  sufficiently  accurate  results.  The  solution 
of  the  coupled  equations  is  achieved  by  iteration  of  approximate  solutions  to  self-consistency. 
The  results  obtained  compare  favorably  to  the  exact  solution  obtained  by  Fuller  (within  5%  for 
size  parameters  less  that  unity). 

The  scattering  and  extinction  cross  sections,  the  differential  scattering  cross  section  and 
the  polarization  ratio  are  calculated  for  agglomerates  consisting  of  17,  52,  and  165  primary 
spheres  for  soot-like  and  silica-like  agglomerates  are  generated  by  a  computer  simulation  of 
Brownian  dynamics  in  flame  growth.  A  comparison  was  made  among  the  Rayleigh-Debye, 
Coupled  Electric  Dipole,  and  Coupled  Electric  and  Magnetic  Dipole  techniques.  It  is  further 
shown  that  the  polarization  ratio  computed  by  the  Coupled  Electric  and  Magnetic  Dipole 
technique  is  sensitive  to  the  primary  sphere  size  independent  of  the  size  of  the  agglomerate 
particle. 
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VARIATIONAL  SCATTERING  TECHNIQUES 

by 

Dr.  Richard  A.  Farrell 
John  Hopkins  Applied  Physics  Laboratory 
Johns  Hopkins  University,  MD 

Variational  techniques  are  a  class  of  solutions  to  electromagnetic  scattering  problems 
which  have  the  advantage  of  quadratic  convergence  to  the  exact  solution.  On  the  other  hand, 
the  technique  is  limited  by  the  accuracy  of  the  initial  guess  of  the  form  of  the  solution.  As  such, 
though  in  principle  the  variational  technique  is  completely  general,  in  practice  it  is  most  useful 
when  the  problem  at  hand  is  an  extrapolation  of  a  problem  whose  solution  is  known  from  other 
techniques.  An  advantage  of  this  approach  is  that  the  variational  technique  can  supply  accurate 
results  at  all  wavelengths.  Of  particular  interest  to  atmospheric  physics  is  the  fact  that  the 
Schwinger  variational  principle  yields  a  formalism  that  is  directly  applicable  to  a  stochastic 
ensemble  of  simple  scatterers  such  as  often  occurs  in  CIAO. 

The  work  discussed  here  extends  the  range  of  validity  of  the  variational  technique  to  long 
wavelengths  via  the  incident  field  approximation  and  to  short  wavelengths  via  the  Kirchoff 
approximation.  It  is  shown  by  comparison  with  measurement  that  this  approach  yields  accurate 
results  for  both  random  and  fixed  agglomerations  of  simple  scatterers,  provided  that  the  trial 
function  takes  proper  account  of  the  refractive  effects  of  dielectrics,  the  boundary  conditions  and 
shadowing  effects  and  the  existence  of  creeping  waves. 

The  paper  shows  conclusively  that  the  variational  technique  can  be  reliably  extended  to 
stochastic  distributions  of  simple  scatterers  and  to  near-field  scattering  problems,  and  that  these 
solutions  are  valid  at  all  wavelength/particle  ratios. 
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COMPOSITE  LIMIT  ANALYSIS  OF  ELECTROMAGNETIC 
CROSS  SECTIONS  OF  ISOMETRIC,  PLANAR  AND  CYLINDRICAL  PARTICLES 

by 

Janon  Embury 

(Presented  by  David  Cohoon) 

U.S.  Army  Edgewood  Reseach 

The  purpose  of  this  research  is  to  study  the  effects  of  particle  size  distributions  and 
particle  shape  distributions  on  extinction.  The  model  developed  is  completely  general,  however 
it  assumes  that  the  concentration  of  the  particulate  is  small  enough  that  multiple  scattering 
effects  within  the  cloud  may  be  ignored. 

Fibers  and  hollow  spherical  shells  are  discussed  in  detail,  as  well  as  a  variety  of 
distribution  functions.  The  model  developed  leads  to  an  understanding  of  the  levels  of  extinction 
that  may  be  expected  for  clouds  of  artificial  aerosols  of  varying  shape,  size,  conductivity  and 
orientation. 
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THE  SCATTERING  OF  ELECTROMAGNETIC  RADIATION  BY  NONCONNECTED, 
HETEROGENEOUS,  BIANISOTROPIC  STRUCTURES 

by 

David  Cohoon 

U.S.  Army  Edgewood  Research 
Development  and  Engineering  Center 
Edgewood,  MD 

Computation  times  and  storage  requirements  limit  the  size  and  complexity  of  problems 
that  can  practically  be  considered.  Volume  integral  equation  techniques  alleviate  this  difficulty 
somewhat.  This  work  describes  a  technique  developed  for  solving  electromagnetic  interaction 
problems  that  can  give  answers  to  machine  precision  by  using  the  concept  of  an  exact  finite  rank 
integral  equation  whose  solution  is  exactly  the  projection  of  the  integral  operator  onto  a  space 
of  approximate  solutions  of  the  original,  infinite  rank,  integral  equations. 

The  integral  equations  discussed  are  valid  for  heterogeneous,  bianisotropic  structures 
which  are  not  necessarily  connected.  It  is  proven  that  this  technique  will  always  converge  to  the 
solution,  but  the  issue  of  how  quickly  it  converges  is  not  settled.  The  technique  is  completely 
general  but  still  new  enough  to  be  untried  on  problems  of  practical  interest.  While  the  method 
is  very  intriguing,  there  is  much  development  work  yet  to  be  done. 
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LIGHT  SCATTERING  BY  INHOMOGENEOUS  SPHERES: 

APPROXIMATE  SOLUTIONS  USING  A  VOLUME  CURRENT  METHOD 

by 

Steven  Hill 
and  Hasan  Saleheen 
New  Mexico  State  University 

It  is  of  interest  to  compute  the  scatter  from  spheres  containing  small  foreign  inclusions 
(See  for  example  the  summary  of  work  presented  by  Ron  Pinnick).  A  technique  called  the 
"volume  current  method"  (VCM)  has  been  developed  to  perform  the  computations.  The  VCM 
assumes  approximate  solutions  for  the  polarization  currents  of  small  inclusions  found  within  the 
primary  sphere,  then  computes  the  resulting  electromagnetic  fields  and  ultimately  the  scattering 
magnitudes  by  integrals  over  the  induced  currents  for  each  inclusion.  The  results  have  been 
compared  with  the  separation  of  variables  code  developed  by  Kirk  Fuller  (see  the  summary  for 
Kirk  Fuller’s  work)  to  establish  a  measure  of  result  reasonableness  and  the  two  methods  are 
shown  to  be  reasonably  consistent.  However,  it  should  be  pointed  out  that  neither  method  has 
been  experimentally  verified. 

Two  comparison  cases  are  considered  between  the  VCM  method  and  the  separation  of 
variables;  nonresonant  spheres,  and  near  resonance  spheres.  For  the  nonresonant  sphere  case, 
an  inclusion  comparable  in  size  to  the  host  is  assumed.  In  this  case  the  VCM  prediction  agrees 
within  5  percent  of  that  predicted  by  the  separation  of  variables  method.  For  the  case  where  the 
inclusion  is  much  less  than  the  host  diameter  and  oriented  within  the  host  for  maximum  scatter 
there  is  almost  perfect  agreement  between  the  two  models. 

When  the  host  sphere  is  near  resonance  in  size,  it  is  important  that  an  effective  average 
refractive  index  be  used  that  accounts  for  both  host  and  inclusions.  The  VCM  and  separation 
of  variables  method  agree  well  in  predicting  the  resonant  peak.  The  scatter  magnitudes  agree 
well. 


It  is  found  that  at  large  scatter  angles  the  differences  between  scatter  from  a  host 
without  inclusions  and  host  with  inclusions  are  large.  This  suggests  that  observations  of  scatter 
magnitude  can  be  used  to  detect  biological  aerosols. 

The  VCM  analysis  conducted  thus  far  assumes  scatter  observed  at  90°  relative  to  the 
incident  beam.  Future  work  should  consider  additional  scatter  angles  and  system  sensitivities. 
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ALGEBRAIC  APPROXIMATIONS  OF  EXTINCTION  EFFICIENCY 
FOR  NON-SPHERICAL  PARTICLES 


by 

Blair  Evans 

Defense  Research  Establishment 
Val  Cartier,  Canada 

The  immediate  objective  of  this  work  is  to  significantly  reduce  the  computational  burden 
in  calculating  the  extinction  from  non-spherical  particles.  The  longer  term  objective  is  to  alleviate 
the  remaining  constraints  in  the  theoretical  consideration  of  non-spherical  aerosols  and 
obscurants  of  either  natural  or  artificial  origin. 

An  analytic  semi-empirical  approximation  for  the  extinction  efficiency  for  randomly 
oriented  spheroids  has  been  developed  and  compared  with  the  extended  boundary  condition  or 
T-matrix  method.  The  basic  approach  is  to  orthogonalize  as  much  as  possible  the  scattering 
physics  into  well-defined  regimes.  For  small  physical  and  optical  sizes,  the  electrostatics 
(Rayleigh)  approximation  is  used.  For  larger  and  very  large  optical  sizes  we  still  use  the 
electrostatics  approximation  but  with  the  optical  constants  transformed  to  include  the  effects  of 
the  magnetic  dipole.  In  the  large  physical  size  regime,  the  analysis  uses  a  diffraction 
(anomalous  diffraction)  component  and  what  can  loosely  be  described  as  edge  effects.  The 
diffraction  component  is  modeled  by  the  anomalous  diffraction  approximation  developed  by  Van 
de  Hulst.  The  edge  effect  (Fock  theory)  component  is  modeled  by  extending  a  technique 
introduced  by  Jones.  This  component  is  further  generalized  to  have  proper  behavior  at  small 
optical  sizes  and  large  indices.  The  gap  between  the  large  and  small  particle  regimes  is  bridged 
by  a  binomial  form  similar  to  the  "generalized  mean." 

The  approximation  has  been  verified  for  complex  refractive  indices,  m  =  n  -  ik,  where  1 
<  n  <  00  and  0  <  k  <  oo  and  aspect  ratios  from  0.2  to  5.  After  extensive  error  analysis  and 
comparisons  with  other  verified  models  such  as  the  T-matrix  and  the  Mie  theory,  it  is  concluded 
that  the  approximation  is  uniformly  valid  over  all  size  parameters  and  aspect  ratios.  Errors 
relative  to  the  T-matrix  range  between  1  and  35  percent.  The  largest  errors  occur  for  k  values 
greater  than  1 ,  and  for  real  refractive  indices  near  1 .  Least  errors,  1  percent  or  less  are  found 
in  comparisons  with  the  Mie  theory. 

All  approximate  computations  for  the  extinction  efficiency  in  this  study  were  produced  at 
a  rate  greater  than  10“  times  faster  than  the  T-matrix  code.  Since  the  T-matrix  computation 
times  scale  at  least  as  the  cube  of  the  optical  size  and  the  analytic  approximation  is  optical  size 
dependent,  larger  size  parameters  or  larger  refractive  indices  lead  to  larger  speed-up  factors. 

Approximation  limitations  for  n<1  and/or  k<0  have  not  been  modeled  since  new  and 
significant  physical  phenomena  arise  (e.g.  total  internal  reflections  and  optical  gain).  Even  n 
modestly  less  than  one  can  cause  problems.  Another  limitation  occurs  for  2  ^  k  <  1 0  and  large 
particles.  In  this  region,  the  absorption  is  not  well  modeled.  If  this  effect  was  properly  accounted 
for,  absorption  and  scattering  efficiencies  could  be  globally  and  readily  obtained  by  using  the 
same  approach. 
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THE  CORPUSCULAR  VIEW  OF  LIGHT 

by 

Mel  Lax 

City  College  of  New  York 

The  corpuscular  view  of  light  promoted  by  Newton  was  discredited  by  Descartes  because 
it  appeared  to  require  an  increase  in  the  velocity  of  light  in  a  dielectric  medium  (using  qualitative 
arguments  based  on  the  direction  of  bending).  The  quantitative  arguments,  as  summarized  by 
Kline  and  Kay  (1965)  and  by  Luneburg  (1966)  seem  to  yield  the  same  conclusion  but  they  are 
based  on  interpreting  an  artificial  parameter  r  as  the  time.  To  reinstate  the  corpuscular  view  of 
light  we  must  find  the  correct  parallelism  between  mechanics  and  light,  not  the  usual  one  chosen 
by  Hamilton. 

Kline,  M.  and  Kay,  I.W.  (1965),  Electromagnetic  Theory  and  Geometrical  Optics,  Wiley- 
Interscience,  Sec.  2.6 

Luneburg,  R.  K.,  (1966)  Mathematical  Theory  of  Optics,  University  of  California  Press. 
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1.3.2  Model  Application  Presentation  Overviews 

Overviews  presented  in  this  section  represent  those  papers  describing  model  applications. 
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CHAFF  MODELING  AND  SIMULATION 

by 

Robert  J.  Puskar,  Ph.D. 

Mission  Research  Corporation 

Electromagnetic  scattering  by  chaff  dipoles  has  been  under  investigation  for  almost  40 
years.  Individual  dipole  scattering  is  well  understood,  scattering  from  clouds  of  dipoles  is  not. 
All  the  dipoles  do  not  scatter  independently  except  for  very  long  times  after  release.  This  means 
linear  addition  of  scatter  from  single  dipoles  is  inappropriate.  Furthermore,  aerodynamic  studies 
show  dipoles  have  preferred  orientations  after  release  into  the  atmosphere.  Thus,  estimating 
scatter  from  randomly  oriented  particles  is  inappropriate.  The  most  accurate  way  to  predict  chaff 
cloud  radar  cross-section  (RCS)  is  to  numerically  solve  the  integral  equation  governing 
scattering.  This  works  for  very  small  chaff  clouds.  Only  small  numbers  of  dipoles  can  be 
handled  in  this  way  because  of  computational  requirements.  However,  the  statistical  properties 
derived  from  this  type  of  analysis  are  valid  for  large  clouds.  Other  scattering  parameters  that 
can  be  modeled  in  this  way  include  frequency  bandwidth,  bistatic  scattering,  and  polarization 
response. 

In  any  chaff  cloud  model,  the  effects  of  dispensing,  dispersal,  blooming,  diffusion,  and 
agglomeration  should  be  included.  No  such  model  exists.  One  of  the  most  detailed  chaff 
aerodynamic  models,  a  6  degree-of-freedom  dipole  trajectory  program  that  predicts  individual 
dipole  trajectories  based  on  particle  aerodynamic  coefficients.  It  does  not  handle  clouds  of 
dipoles.  To  a  large  extent,  this  aspect  of  chaff  cloud  modeling  relies  on  empirically  derived  data. 

There  are  four  basic  types  of  chaff  models  in  use  in  the  electronic  warfare  community. 
The  first  type  is  a  chaff  signature  simulation  and  is  used  to  analyze  and  study  chaff  radar 
signatures  in  fine  detail.  The  second  is  a  chaff/radar  effectiveness  model  and  is  used  primarily 
to  evaluate  the  effects  of  chaff  on  radar  system  performance.  The  third  and  fourth  types  are 
used  to  study  the  aerodynamic  and  electromagnetic  scattering  characteristics  of  chaff 
respectively.  Nine  different  models  used  in  these  categories  are  described  in  this  work. 

There  are  three  primary  modes  in  which  chaff  is  now  used:  saturation,  self-protection,  and 
confusion/deception.  The  way  chaff  is  used  can  be  extremely  important  in  determining 
effectiveness,  especially  in  the  self-protection  mode  which  is  now  the  most  common  use  for 
chaff.  Saturation  chaff  refers  to  the  practice  of  producing  a  huge  corridor  of  chaff  to  mask  the 
penetration  of  other  aircraft.  This  practice  is  not  likely  to  be  used  in  the  future  because  of  tactics 
now  used  for  penetration  by  the  USAF.  Confusion/deception  refers  to  the  generation  of  false 
radar  targets.  The  effectiveness  of  this  mode  is  very  difficult  to  evaluate.  There  no  longer 
seems  to  be  significant  effort  by  the  USAF  to  use  chaff  in  a  confusion/deception  mode.  In  an 
Army  battlefield  situation  it  could  have  major  effects. 


30 


BACKSCATTERING  BY  NONSPHERICAL  HYDROMETERS  AS  CALCULATED  BY  THE 

COUPLED  DIPOLE  METHOD:  AN  APPLICATION  IN  RADAR  METEOROLOGY 

by 

Clifton  E.  Dungey 

Department  of  Engineering  Physics 
Air  Force  Institute  of  Technology 
Wright-Patterson  AFB  OH  45433-7765 

and 

Craig  F.  Bohren 
Department  of  Meteorology 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania  16802 

The  ability  of  the  coupled  dipole  method  (also  called  the  "Discrete  Dipole  Approximation" 
by  Flatau)  to  calculate  backscattering  by  nonspherical  particles  is  demonstrated.  For  spheres 
smaller  than  the  size  which  produces  the  first  backscattering  minimum  as  a  function  of  size 
parameter,  small  shape  variations  in  the  dipole  array  do  not  affect  the  accuracy  of  the 
calculation.  Conversely,  when  spheres  are  larger  than  the  size  which  produces  the  first 
minimum,  backscattering  calculations  lose  accuracy.  For  particles  with  an  aspect  ratio  of  at  least 
5:1  and  small  size  parameter,  the  backscattering  calculation  is  relatively  insensitive  to  the  exact 
cross-sectional  shape  of  the  particle,  i.e.  hexagonal  column,  cylinder,  etc.  Backscattering 
calculations  for  ice  spheres,  hexagonal  columns,  and  plates  are  presented.  Polarization  effects 
are  analyzed  using  principles  associated  with  the  Mueller  matrix. 

We  find  identification  of  ice  crystal  size  distribution  depends  on  several  factors.  First,  the 
size  distribution  is  assumed  continuous.  A  discontinuous  size  distribution  produces  gaps  in  the 
Doppler  shift  spectrum  used  for  backscatter  analysis  that  could  be  mistaken  for  a  backscattering 
minimum.  Second,  for  our  theory  to  be  applicable,  ice  crystals  must  not  be  aggregated. 
Backscatter  computations  have  only  been  made  for  plates  and  columns;  scattering  by 
aggregates  may  obscure  the  desired  backscattering  minima.  Third,  the  optical  depth  of  the  ice 
cloud  must  not  be  too  large.  Othenwise,  the  back  scattered  signal  will  be  a  result  of  multiple 
scattering  events.  Our  procedure  for  identifying  ice  crystals  is  based  on  single  scattering. 

Our  results  show  it  is  feasible  to  estimate  size  distribution  for  hexagonal  columns  and 
plates  using  millimeter  wave  radar.  In  general,  the  size  of  the  ice  crystals  dictates  the 
wavelength  best  suited  for  obtaining  the  first  backscattering  minimum.  However,  at  shorter 
wavelengths  attenuation  by  atmospheric  water  vapor  may  limit  system  performance. 
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LIGHT  SCATTERING  BY  AGGREGATES  AND 
INHOMOGENEOUS  SPHERES 

by 

Kirk  A.  Fuller 

Department  of  Atmospheric  Science 
Colorado  State  University 

Electromagnetic  scattering  by  agglomerations  of  spheres  or  by  spherical  hosts  possessing 
such  agglomerates  or  eccentrically  located  spheres  can  be  solved  within  a  single  theoretical 
framework.  The  theory  uses  a  repeated  application  of  the  same  principles  employed  in  the 
Lorenz-Mie  treatment  of  light  scattering  by  single  homogeneous  spheres.  The  approach  is  to 
first  transform  a  plane  wave  incident  on  a  spherical  particle  into  vector  spherical  harmonic  waves 
that  can  be  matched  directly  to  the  spherical  geometry.  The  Mie  coefficients  for  the  expansion 
can  then  be  interpreted  as  being  analogous  to  Fresnel  reflection  and  transmission  coefficients 
for  a  plane  wave  incident  on  a  planar  slab  of  material.  To  extend  the  theory  to  an  ensemble  of 
spheres  of  arbitrary  number  and  size,  the  scattered  fields  from  each  sphere  are  added 
coherently  by  using  the  addition  theorem  for  vector  spherical  harmonics  and  applying  the  proper 
boundary  conditions  for  each  particle.  This  approach  requires  the  solution  of  a  set  of  coupled 
equations  and  represents  the  practical  computational  limit  on  the  number  of  particles  in  an 
aggregate  which  can  be  evaluated. 

The  procedure  has  been  applied  to  the  computation  of  the  extinction  cross-section  of 
small  carbon  particles  on  the  surface  of  sulfate  droplets.  The  extinction  cross-section  varies  with 
the  carbon  particle  position  on  the  droplet  surface  relative  to  the  plane  of  incidence  and  the 
incident  polarization.  When  the  particle  in  or  near  the  focal  volume  of  the  drop  scattering  is 
enhanced,  when  the  particle  is  removed  from  the  focal  volume  area  scatter  is  reduced.  For  an 
average  over  all  orientations  and  polarizations,  the  absorption  characteristics  of  the  carbon 
particles  do  not  change  much  relative  to  the  free  particles.  Placing  the  carbon  grain  inside  a 
droplet  shows  that  the  absorption  cross-section  decreases  as  the  particle  moves  from  the  center 
of  the  drop  to  the  surface  when  averaged  over  all  polarizations  and  orientations. 

When  the  procedure  is  applied  to  a  linear  chain  of  spheres  and  an  average  is  performed 
over  all  orientations  and  incident  polarizations,  the  ratio  of  aggregate  to  single  equivalent  sphere 
cross-sections  increases  by  about  50%  for  carbon  particles  and  about  300%  if  the  material  is 
aluminum. 

The  solution  procedure  developed  for  computing  cross-sections  for  aggregates  of  spheres 
is  very  powerful  and  accurate.  For  example,  it  yields  a  solution  that  is  much  simpler  than 
previously  computed  when  applied  to  the  coated  sphere  problem.  The  analysis  yields 
expressions  in  terms  of  the  Mie  coefficients  solution  for  a  single  sphere  that  are  analogous  to 
the  reflection  and  transmission  coefficients  for  a  thin  film.  The  procedure  for  solving  the 
extinction  from  aggregates  provides  a  method  for  testing  the  accuracy  of  approximations  used 
for  solving  the  scatter  from  irregularly  shaped  particles. 
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COMPUTATION  OF  RADIATIVE  CROSS  SECTIONS 
OF  SPHERE  CLUSTERS 

by 

Daniel  MacKowski 
Auburn  University 

The  motivation  for  this  work  is  the  modeling  of  absorption,  scattering,  and  'snction  of 
super  spherical  agglomerates  such  as  carbonaceous  soot  arising  during  combustion  processes. 
The  analytical  formulation  of  the  solution  to  Maxwell’s  equations  for  this  problem  is  similar  to  that 
of  Kirk  Fuller  where  a  superposition  of  vector  spherical  harmonic  expansions  about  each  of  the 
spheres  is  used.  The  analysis  shows  that  simple  recurrence  relationships  can  be  used  to  obtain 
the  vector  spherical  harmonic  addition  coefficients.  It  is  found  for  this  problem  that  it  is  simpler 
to  compute  the  scattering  cross-section  from  the  difference  of  the  extinction  and  absorption 
cross-sections,  than  is  normally  done  by  computing  the  scattering  and  extinction  cross-sections 
first.  Random  particle  orientations  can  be  addressed  simply  by  using  a  T  matrix  (not  the  T- 
matrix  for  extinction  computations)  that  is  developed  in  the  problem  formulation.  Comparisons 
of  analytical  predictions  with  previously  acquired  data  show  agreement  within  about  30  percent 
for  the  extinction  efficiency. 

For  visible  and  infrared  wavelengths,  and  size  parameters  much  less  than  1, 
computational  results  show  that  more  than  just  the  dipole  terms  are  needed  for  the  computation. 
For  example,  computations  for  a  binary  cluster  with  size  parameter  of  0.01,  shows  that  as  the 
complex  index  of  refraction  increases  so  does  the  number  of  terms  required  to  specify  the 
absorption  efficiency.  However,  using  the  electrostatic  approximation  (Rayleigh  approximation) 
for  the  same  case  shows  that  the  absorption  efficiency  is  constant  for  complex  index  of  refraction 
when  the  size  parameter  is  less  than  about  0.1 ,  which  is  as  expected.  A  similar  problem  is  found 
for  particles  consisting  of  large  numbers  of  small  size  parameter  spheres  in  an  aggregate.  It  is 
believed  that  the  differences  in  the  analytical  results  arise  because  a  pair  of  dipoles  act  together 
more  like  a  spheroid,  than  a  simple  dipole  pair. 

An  approximation  for  the  extinction  by  spherical  aggregates  has  been  developed  by 
considering  interactions  only  between  pairs  of  dipoles.  The  approximation  gives  a  closed  form 
solution  that  is  a  slight  under  prediction  of  the  absorption  cross-section  computed  for  multiple 
particle  interactions  beyond  pairs  only.  An  appropriate  characteristic  size  of  the  aggregate  is 
found  to  be  the  radius  of  gyration  divided  by  the  number  of  fundamental  particles  in  the 
aggregate  raised  to  the  1/3  power. 
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ORIENTATION  EFFECTS  OF  FALLING 
NON-SPHERICAL  PARTICLES 

by 

James  Klett 
PAR  Associates 

Prediction  of  falling  non-spherical  particle  orientation  is  addressed  by  a  heuristic  model 
that  assumes  the  particles  are  subject  to  isotropic  turbulence  within  or  below  the  inertial 
subrange,  i.e.  the  Kolmogorov  spectrum  of  eddies,  depending  on  particle  size.  The  RMS  tilt 
angle  of  a  small  eccentricity  spheroidal  particle  is  determined  by  Langevin-type  averaging  over 
its  equation  of  motion.  The  calculation  takes  into  account  the  first  order  restoring  torque  arising 
when  the  stable  fall  mode  is  perturbed  by  either  thermal  or  turbulent  fluctuations.  An 
approximate  expression  for  the  variance  of  an  assumed  Gaussian  orientation  distribution  for 
small  tilt  angles  and  small  Reynolds  numbers  is  obtained  by  invoking  dimensional  constraints 
concerning  the  nature  of  the  main  flow  and  turbulent  stresses,  and  by  assuming  turbulent  and 
thermal  fluctuations  are  uncorrelated.  The  expression  is  then  generalized  to  provide  a  semi- 
quantitative,  nearly  Gaussian  orientation  distribution  for  arbitrary  tilt  angles,  particle  aspect  ratios, 
Reynolds  numbers,  and  particle  sizes  relative  to  the  Kolmogorov  microscale  length.  The  analysis 
applies  to  particles  that  can  be  modeled  as  spheroids,  disks,  and  cylinders,  as  well  as  ice  crystal 
hexagonal  plates  and  columns. 

The  principal  analytical  result  is  an  orientation  distribution  of  the  form  p(0)  =  2ifi)  exp(- 
/ffsin^(0)),  where  Z  is  the  normalization  constant,  and  where  p  is 

e  is  the  tilt  angle  measured  from  the  stable  position,  y  is  the  particle  aspect  ratio,  m  is  the 
effective  mass  term  representing  the  mass  of  air  in  a  volume  about  54  percent  of  the  particle 
volume,  Npe  is  the  particle  Reynolds  number,  k  is  a  number  of  0(1),  is  the  particle  terminal 
speed,  k  is  Boltzmann’s  constant,  T  is  absolute  temperature,  is  the  Reynolds  number  of 
turbulent  eddies  of  scale  size  equal  to  the  particle  maximuni  dimension,  and  v,  is  the 
corresponding  magnitude  of  the  turbulent  eddy  fluctuating  velocity.  The  kT  term  in  the 
denominator  of  represents  the  disorienting  effects  of  Brownian  motion,  and  the  limiting  form 
of  the  distribution  for  the  case  of  vanishing  turbulence  coincides  with  the  distribution  due  to 
thermal  effects  obtained  by  Fraser  (JOSA,  69, 1 1 1 2  (1 979)).  According  to  the  model,  turbulence 
can  significantly  affect  the  orientational  order  of  some  aerosols.  For  example,  for  the  case  of 
carbon  and  glass  fibers  of  3mm  length  and  density  2.5g/cm®,  falling  in  mildly  turbulent  air 
(characterized  by  a  turbulent  energy  dissipation  rate  of  10  cm^/sec^ ),  a  transition  from  orderly 
fall  to  random  tumbling  is  predicted  if  the  diameters  of  the  fibers  are  reduced  from  1 0//m  to  5^m. 
This  result  is  consistent  with  field  observations  (W.  Michael  Farmer,  private  communication, 
1991). 
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SIMPLE  MODEL  AND  MEASUREMENTS  FOR  MODEL  VALIDATION  IN  THE  MM/CM 

WAVELENGTH  REGION 

by 

Ira  Goldberg 

Rockwell  International  Space  Center 

A  simple  semi-empirical  approach  is  developed  to  evaluate  the  attenuation  and  scattering 
in  the  mm  and  cm  wavelength  regions  of  particulate  dispersed  in  air.  The  approach  involves  an 
improved  Maxwell-Gamett  type  model  to  calculate  the  attenuation  as  a  function  of  frequency. 
The  ensemble  of  particles  is  treated  as  an  artificial  dielectric  in  which  the  aspect  ratio  and  the 
conductivity  are  the  input  variables.  These  calculations  are  compared  to  data  measured  on  a 
suspension  of  ellipsoids.  Several  materials  have  been  studied:  carbon  cylinders,  nickel  coated 
carbon,  aluminum  coated  glass  and  iron  coated  glass.  The  best  fit  to  the  experimental  data  is 
obtained  by  using  the  measured  effective  conductivity  of  the  fibers  rather  than  the  bulk 
conductivity.  The  effective  conductivity  is  generally  about  10%  lower  than  the  bulk  conductivity. 

Validation  is  done  by  measuring  the  absorption  and  scattering  by  a  fluid  suspension  of 
filaments  uniformly  dispersed  in  a  fixed  volume  tank.  Measurements  are  taken  in  the  range  from 
2-20  GHz  and  from  26-140  GHz.  The  measurement  technique  allows  the  concentration 
dependence  to  be  accurately  determined  and  effects  of  random  versus  aligned  orientation  to  be 
considered.  Experience  has  shown  that  the  experimental  setup  is  reproducible. 

The  main  conclusion  drawn  from  this  work  is  that  the  modified  Maxwell-Gamett  model 
provides  a  quick  and  accurate  computational  tool  whose  predictions  agree  well  with  the 
measured  performance. 


35 


MM  WAVE  POLARIZATION  FEATURES 
FOR  AUTOMATIC  TARGET  RECOGNITION  AND  DETECTION 

by 

William  Cameron 
The  Boeing  Co 
Seattle,  WA 

MM  wave  target  signatures  are  recognized  by  spatially  resolving  the  constituent 
scatterers.  The  maximum  degree  of  discrimination  is  achieved  by  using  polarization  effects  to 
identify  the  individual  scatterer  types  and  orientations.  These  scatterer  types  and  orientations 
are  extracted  from  the  polarization  scattering  matrix  by  transforming  it  into  a  vector  for 
subsequent  analysis.  It  is  shown  that  for  symmetric  scatterers,  the  responses  can  be 
represented  as  points  contained  in  the  unit  disc  of  the  complex  plane. 

These  methods  are  applied  to  measured  signatures  of  simple  polarimetric  scatterers  to 
determine  the  scatterer  type  and  orientation.  Several  types  of  examples  are  discussed  in  detail. 
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1.3.3  Experimental  Evaluation  and  Analyses  Overviews 

Overviews  presented  in  this  section  represent  those  papers  describing  experimental  evaluation 
and  analyses  of  models. 
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THE  SCATTERING  OF  INFRARED  RADIATION 
BY  IRREGULAR  PARTICLES  ON  SURFACES 

by 

Lowell  D.  Lamb.  Julie  D.  Lorentzen,  and  Donald  R.  Huffman 
Department  of  Physics,  University  of  Arizona 
Tucson,  Arizona 

Mueller  scattering-matrix  elements  were  measured  for  50  irregular  SiC  particles  on  a 
NaCI  substrate.  The  measurements  were  made  with  a  polarization-modulated,  angle-scanning, 
infrared  scattering  instrument  that  has  a  tunable  wavelength  COj  laser  as  its  source  of 
illumination.  The  particles  were  illuminated  individually  at  two  wavelengths  (9.25  and  10.66//m), 
and  a  small  angular  region  of  the  scattered  light  was  collected  by  a  Cassegrain  telescope  held 
at  a  fixed  angle.  Because  the  electromagnetic  behavior  of  SiC  changes  dramatically  as  the 
illumination  wavelength  is  varied  from  "insulating"  at  9.25//m  to  "metallic"  at  10.66//m,  it  was 
possible  to  measure  the  scattering  of  pairs  of  systems  which,  except  for  particle  composition, 
were  identical  in  every  way.  Hence,  the  need  to  physically  remove  one  particle  and  replace  it 
with  another  and  the  many  potential  sources  of  experimental  error  which  would  accompany 
replacement  of  the  particle  were  eliminated.  Comparison  of  a  statistical  analysis  of  the  data  to 
a  sum-of-fields  model  suggests  that  measurements  of  the  polarization  states  of  light  scattered 
from  irregular  particles  on  surfaces  could  be  used  potentially  in  practical  applications  to 
discriminate  between  compositional  classes  of  particles. 
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LIGHT  SCATTERING  FROM  ASSEMBLIES  OF  SCATTERS: 

APPLICATIONS  TO  THE  CORNEA  OF  THE  EYE 

by 

D.E.  Freund.  R.L.  McCally  and  R.A.  Farrell 
The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  MD  20723 

When  making  comparisons  between  theory  and  experiment,  care  must  be  taken  to  assure 
that  the  quantity  explicitly  computed  corresponds  to  what  is  measured,  and  that  the  analysis 
model  does  not  neglect  any  of  the  relevant  physical  processes  governing  the  experimental  data. 
Two  examples  based  on  our  research  using  scattered  light  to  probe  the  ultrastructure  of  the 
cornea  are  presented  to  illustrate  the  importance  of  these  considerations.  We  show  how  it  is  not 
always  easy  to  recognize  when  one  or  the  other  is  being  violated. 

The  first  example  shows  the  importance  of  computing  quantities  that  correspond  to  what 
is  measured.  In  particular  it  shows  that  calculations  performed  over  an  extremely  small 
subregion  give  rise  to  a  diffraction  term  that  depends  only  on  the  size  and  shape  of  the 
subregion  and  not  on  the  scatterers  contained  within.  The  computed  results  (visible 
transmittance  through  the  cornea)  disagreed  with  experimental  data  that  probed  a  macroscopic 
region  in  which  only  the  scattering  from  the  assembly  of  scatterers  was  measured.  However, 
by  recomputing  the  scattered  intensity  to  more  nearly  match  the  experimental  conditions  under 
which  the  data  were  acquired  by  eliminating  the  diffraction  term,  we  obtained  excellent 
agreement  between  theory  and  experiment. 

The  second  example  shows  the  importance  of  assuring  the  idealized  model  contains  all 
the  relevant  physical  processes.  In  this  example  the  model  gives  good  results  for  total  scattering 
cross-section  but  poor  results  for  the  angular  scattering  cross-section.  This  example  is 
particularly  striking  because  the  total  scattering  cross-section  was  obtained  by  first  computing 
the  angular  scattering  cross-section  and  then  numerically  integrating  to  obtain  the  total  scattering 
cross-section.  We  found  that  only  when  the  model  was  refined  to  include  highly  significant 
orientation  effects  of  the  constituent  scatterers  did  theoretical  predictions  agree  with  the 
measurements. 
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PARAMETRIC  SENSITIVITY  STUDIES  OF  A  VARIATIONAL  MODEL 
FOR  MILLIMETER  WAVE  OBSCURANTS 

by 

W.  Michael  Farmer 
Bruce  Kennedy 

The  Physical  Sciences  Laboratory 
New  Mexico  State  University 

A  sensitivity  analysis  was  performed  on  the  commonly  used  Pedersen-Waterman 
"extended  quasistatic"  model  to  predict  the  extinction  cross-section  for  high  aspect  ratio 
cylinders.  The  analysis  was  applied  to  Dixon  1102  graphite  and  aluminum  millimeter  wave 
obscurants.  The  objectives  of  the  study  were  to  determine  the  impact  of  uncertainty  in  input  data 
on  model  prediction  capability  and  to  identify  the  most  significant  experimental  parameters  that 
affect  model  prediction.  Study  results  show  that  model  uncertainty  is  strongly  dependent  on 
particle  material,  and  is  relatively  insensitive  to  variations  in  particle  length.  For  the  particle 
diameters  and  lengths  considered  in  this  study,  10%  relative  uncertainty  in  dielectric  constants 
and  length  for  graphite  was  found  to  yield  higher  relative  uncertainty  in  extinction  cross-section 
predictions  at  frequencies  below  about  60  GHz  than  aluminum.  Above  60  GHz,  aluminum 
extinction  cross-section  predictions  yielded  higher  relative  uncertainties  than  graphite. 
Uncertainty  in  model  predictions  are  near  the  uncertainties  in  particle  length,  all  other  tactors 
having  no  relative  uncertainty.  For  aluminum,  relative  uncertainty  in  model  predictions  is 
essentially  due  to  the  imaginary  component  of  the  dielectric  constant  for  frequencies  below 
40GHz.  For  graphite,  the  real  and  imaginary  components  make  nearly  equal  contributions  for 
frequencies  as  high  as  10GHz,  and  then  become  dominated  by  the  imaginary  component  of  the 
dielectric  constant.  The  relative  uncertainty  in  transmission  as  a  result  of  uncertainties  in 
extinction  cross-section  predictions  was  found  to  increase  as  transmission  decreases  and  to 
increase  as  frequency  decreases. 

Study  results  clearly  demonstrate  that  model  predictions  are  generally  only  as  accurate 
as  input  values.  The  analysis  shows  that  model  prediction  uncertainty  typically  varies  between 
1.5  and  2.5  times  that  of  the  input  parameters.  However,  for  frequencies  in  the  60-100GHz 
range,  relative  uncertainty  in  extinction  cross-sections  for  aluminum  due  to  relative  uncertainty 
in  dielectric  constant  are  significantly  less  than  the  relative  uncertainty  of  the  input  values. 

For  reliable  model  predictions  it  is  important  that  relative  uncertainty  in  the  prediction  be 
estimated  as  a  function  of  the  relative  uncertainty  of  model  input  parameters.  The  model  may 
be  a  perfect  representation  of  physical  reality  and  still  fare  poorly  in  its  predictive  capability 
simply  because  of  the  uncertainty  in  model  inputs  and  because  of  the  model’s  sensitivity  to  those 
inputs. 
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SUPPRESSION  OF  SCATTERING  RESONANCES  IN  INHOMOGENEOUS 

MICRODROPLETS 

by 

Ron  Pinnick 

U.S.  Army  Research  Laboratory 
Battlefield  Environment  Directorate 
White  Sands  Missile  Range,  NM  88002-5501 

This  work  considers  the  prediction  and  experimental  measurement  of  the  scatter  from 
liquid  droplets  containing  much  smaller  randomly  distributed  spherical  inclusions.  Particles  of 
this  type  include  atmospheric  hydrometers  that  have  scavenged  small  solid  aerosols,  e.g.  sulfate 
based  droplets  containing  carbon  particles,  dust  particles  with  individual  mineral  inclusions,  and 
liquid  droplets  bearing  biological  materials.  It  is  this  latter  case  which  is  the  primary  objective 
for  analysis  in  this  study. 

For  this  study,  glycerol  droplets  2-10  //m  in  diameter  are  levitated  and  suspended  in  an 
electric  field  trap  and  illuminated  with  a  laser.  Scatter  from  the  droplet  is  observed  at  90° 
relative  to  the  incident  beam.  The  droplet  diameter  slowly  changes  through  evaporation  over 
observational  periods  of  minutes  to  hours.  The  droplet  diameter  is  precisely  sized,  with  an 
uncertainty  of  less  than  ±10nm,  using  scattering  resonances  predicted  by  Mie  theory.  Theory 
and  experiment  are  in  sufficient  agreement  that  over  80  observed  resonance  cycles  can  be 
matched  with  theoretical  predictions. 

When  the  glycerol  droplet  is  impregnated  with  10^  to  10®  latex  spheres  having  diameters 
of  the  order  of  60nm,  the  resonance  amplitude  decreases  and  is  broadened.  These  effects  are 
observed  to  increase  in  magnitude  as  the  size  of  the  added  latex  inclusions  increases.  The 
signal  is  also  observed  to  fluctuate  in  time  with  the  presence  of  latex  inclusions,  whereas  for  the 
pure  glycerol  drop  the  signal  does  not  have  the  temporal  fluctuation  component. 

In  an  attempt  to  explain  the  experimental  observations,  the  computer  code  developed  for 
predicting  scatter  by  spherical  aggregates  (see  the  summary  of  Kirk  Fuller’s  work)  was  used  to 
compute  scatter  effects  by  a  single  small  inclusion  as  it  drifted  through  high  intensity  illumination 
regions  inside  the  glycerol  droplet.  The  scatter  predicted  for  this  case  is  insufficient  to  explain 
the  experimental  results.  To  obtain  sufficient  scatter  magnitude  would  require  an  agglomeration 
of  about  11  latex  spheres.  This  is  not  believed  to  be  likely  to  occur.  A  second  possible 
explanation  considers  coherent  scatter  by  pairs  of  latex  particles  that  produce  homodyne  Doppler 
shift  interference  that  temporally  modulates  the  signal.  Analysis  of  this  effect  is  in  process. 
Regardless  of  the  correct  physical  explanation,  the  observed  phenomena  suggest  the  possibility 
that  resonance  scatter  has  the  potential  for  distinguishing  particles  with  and  without  inclusions. 
This  is  a  technique  worth  consideration  for  the  detection  of  biological  aerosols. 
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1.  General  Millimeter  Wave  Obscurant  Bibliography 

The  following  bibliography  represents  reports  and  refereed  journal  publications  that  have  been 
i>alected  to  provide  a  reasonable  overview  of  the  state  of  data  sources  and  data  bases  for  data 
with  which  to  evaluate  electromagnetic  wave  scattering  and  propagation  models  not  generally 
published  in  the  open  literature.  The  bibliography  is  divided  into  sections  that  parallel  those 
discussed  in  section  5.  Many  of  the  technical  report  references  listed  for  millimeter  wave  and 
microwave  chaff  were  extracted  from  the  bibliography  assembled  by  Gary  Grider  at  Wright- 
Patterson  Air  Force  Base. 
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2.  ISST  Data  Base  Summary 

This  appendix  contains  a  brief  summary  of  scattering  measurements  made  at  the  Institute  of 
Space  Science  and  Technology,  Gainsville,  FL.  As  provided  by  Dr.  Ru  Wang. 

TABLE  B-1.  Selected  Microwave  Scattering  Data  for  Rough  Particles  and  the  Target 
Parameters  - 


ID 

# 

Target 

Shape 

do(cm) 

Outer 

Diam. 

E83 

Longer 

Length 

Ucm) 

Shorter 

Length 

Xv 

Volume 

Equiv. 

Size 

Refractive 

Index 

Extinc¬ 

tion 

P-Q 

Data 

Angular 

Scatt 

Data 

PO  RO 

1 

A 

3.782 

1.891 

3.777 

1.885 

3.650 

1.610-/0.004 

14 

34 

14 

2 

D 

5.692 

2.856 

5.692 

2.843 

5.499 

1.610-/0.004 

14 

34 

14 

3 

7.607 

3.809 

7.620 

3.793 

7.346 

1.610-/0.004 

14 

34 

14 

4 

A 

8.905 

4.482 

8.956 

4.482 

8.626 

1.610-/0.004 

14 

34 

14 

5 

A 

3.439 

2.144 

3.490 

2.228 

3.592 

1.256-/0.003 

14 

34 

14 

6 

A 

5.354 

2.702 

5.377 

2.686 

5.184 

1.388-/0.005 

14 

34 

14,V 

B 

B 

1.529 

4.580 

3.053 

4.257 

1.366-/0.005 

14 

34 

14,4 

8 

B 

2.187 

6.538 

4.314 

6.069 

1.367-/0.005 

14 

34 

14.V 

9 

B 

2.754 

8.161 

5.582 

7.694 

1.360-/0.005 

14 

34 

14 

10 

B 

4.622 

13.854 

9.260 

12.882 

1.353-/0.005 

14 

34 

14 

Legend: 

(1)  Target  Shape: 


Shape  A 

—  W 


4> 


— 1, 


tog- 


Shape  B 

— W  — 


(2)  PO:  Data  for  a  Preferentially  Oriented  particle.  RO:  Averaged  data  over  Radom  Orientations. 

(3)  References  where  the  marked  data  are  cited  or  described: 

♦:  Wang,  R.  T.,  1984  in  Proc  (1)  ♦:  Wang,  R.  T.,  1984,  in  Proc  (2) 

v:  Wang,  R.  T.,  1992,  "Scattering..."  ♦:  Wang,  R.  T.,  1992,  "Extinction..." 
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Table  B-2.  Selected  Microwave  Scattering  Data  for  Finite  Cylinders  and  the  Target 
Parameters 


ID 

# 

a 

radius 

(cm) 

/ 

length 

(cm) 

p^l/(2a) 

aspect 

ratio 

ka-2n  a/A 
cylinder 
size 

parameter 

Xv 

volume- 

equival. 

size 

param. 

refractive 

index 

Extinc¬ 

tion 

P-Q 

Data 

An( 

Sc 

D 

PO 

3ular 

:att. 

ata 

RO 

Scatt. 
Intensity 
vs.  Azim. 
Sweep 

1 

.849 

3.400 

2.002 

1.676 

2.418 

1.610-/.  004 

1 

2 

.955 

3.828 

2.004 

1.885 

2.720 

II 

1 

3 

1.058 

4.240 

2.003 

2.089 

3.014 

M 

1 

B 

1.165 

4.670 

2.004 

2.299 

3.318 

It 

1 

1.270 

5.100 

2.008 

2.506 

3.620 

II 

1 

3 

1 

6r 

6 

1.376 

5.512 

2.002 

2.717 

3.920 

II 

1 

7 

1.482 

5.940 

2.003 

2.926 

4.222 

II 

1 

8 

1.590 

6.365 

2.002 

3.137 

4.526 

II 

1 

9 

1.799 

7.206 

2.003 

3.551 

5.123 

II 

1 

10 

2.012 

8.050 

2.001 

3.970 

bill 

II 

1 

11 

2.118 

8.485 

2.004 

4.179 

6.031 

II 

1 

12 

2.328 

9.328 

2.004 

4.594 

6.630 

II 

1 

13 

2.540 

10.185 

2.004 

5.014 

7.237 

II 

1 

14 

1.304 

5.014 

1.922 

2.574 

3.663 

1.322-/.  005 

1 

15 

1.302 

5.032 

1.932 

2.570 

3.671 

1.4 16-/.  005 

1 

16 

2.490 

10.056 

2.019 

4.914 

7.112 

1.321-/.005 

1 

17 

2.512 

10.046 

2.000 

4.958 

7.118 

1.442-/.005 

1 

■■ 

18 

.785 

6.282 

4.001 

1.549 

2.816 

1.610-/.004 

1*,4 

3 

19 

.964 

7.71'2 

4.000 

1.903 

3.457 

II 

1* 

3 

1 

20 

1.204 

9.636 

4.002 

2.376 

4.320 

1*,4 

3 

1 

21 

1.457 

11.656 

4.000 

2.876 

5.225 

II 

14,« 

3 

1 

22 

1.157 

9.256 

4.000 

2.284 

4.150 

1.335-/.  005 

1 

3 

1 

6r,+ 

23 

1.620 

12.970 

4.002 

3.198 

5.812 

1.390-/.  005 

1 

■I 

■m 

■■■ 

24 

.635 

10.110 

7.961 

1.253 

2.816 

1.610-/.  004 

1 

BE 

25 

1.144 

18.334 

8.009 

2.259 

5.174 

1.333-/.  005 

1 

26 

.744 

18.260 

12.26 

1.469 

3.878 

1.327-/.  005 

1 

34 

1 

6r,* 

Legend: 

(1)  PO:  Data  for  a  preferentially  oriented  cylinder. 

RO:  Averaged  data  over  randon  orientaions. 

(2)  References  where  the  marked  data  are  cited  or  described: 

:  Schuerman,  D.  W.  ,  R.  T.  Wang  ,  1981  ♦:  Wang,  R.  T.  and  D.  W.  Schuerman,  1981 
♦:  P.  e.  Geller,  1985  Wang,  R.  T.,  and  Y.  L.  Xu,  1985 

4:  Wang,  W.  X.,  and  R.  T.  Wang  +:  J.  I.Hage,  1991 
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Table  B-3.  Selected  Microwave  Scattering  Data  for  Circular  Disks  &  Square  Plate  and  the 
Target  Parameters 


ID 

# 

a(cm) 
radius 
or  edge 
length 

thick¬ 

ness 

P= 

tl{2a) 

aspect 

ratio 

ka= 

2nalA 

X, 

volume- 
equival. 
size  param 

m=/T7-/m" 

refractive 

index 

Extinc¬ 

tion 

P-Q 

Data 

Angular 

Scattering 

Data 

PO  RO 

Scatt. 

Intensity 

versus 

Azimuth 

1 

1.8978 

0.4711 

0.1241 

3.746 

2.139 

1.610-i.004 

1 

■1 

2 

1.8862 

0.9484 

0.2514 

3.723 

2.690 

1.380-1005 

1 

3 

1.924 

0.962 

0.250 

3.797 

2.738 

1.610-1.004 

1*,4 

3 

1 

D 

2.361 

1.180 

0.250 

.  4.660 

3.360 

1.610-1.004 

1*.4 

3 

1 

2.950 

1.476 

0.250 

5.822 

5.123 

1.610-i.004 

1*,4 

3 

1 

6 

3.569 

1.784 

0.250 

7.044 

6.198 

1.610-i.004 

1^,4 

3 

1 

7 

2.497 

2.537 

0.500 

4.928 

4.501 

1.318-1.005 

1 

8 

(Square  Plate):  4.7612cm  x  4.7612cm 

X  0.9492cm 

3.405 

1.366-1.005 

ir 

Legend: 

(1)  PO:  Data  for  a  preferentially  oriented  disk. 

RO:  Averaged  data  over  random  orientations. 

(2)  References  where  the  marked  data  are  cited  or  described: 

:  Schuerman,  D.  W.,  R.  T.  Wang  ,  1981 
#:  Wang,  R.  T.,  and  D.  W.  Schuerman;  Wang,  R.  T.,  1982  "Extinction..." 

♦  :  Wang,  W.  X.,  and  R.  T.  Wang,  1982  "Corrections..." ;  Wang,  W.  X.,  and  R.  T.  Wang,  1987  "Theoretical..." 
y:  Hage,  J.  I.,  1991 
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Appendix  C 


Results  From  CiAO  II 
Modeling  Exercise 
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1.  Purpose  and  Scope 

This  appendix  summarizes  the  results  of  model  computations  made  by  CIAO  II  participants 
using  the  same  aerosol  sizes  and  materials  complex  index  of  refraction.  The  objective  of  this 
comparison  was  to  determine  how  well  different  models  computing  the  same  electromagnetic 
wave  scattering  parameters  agreed  and  to  identify  those  areas  where  caution  needs  to  be 
exercised  in  the  use  of  these  models. 

Models  developed  for  millimeter  wave  scattering  and  absorption  by  fibers  and  for  infrared 
extinction  by  disks  were  exercised  prior  to  the  CIAO  II  workshop.  Results  were  tabulated  and 
compared  for  illumination  by  three  millimeter  wave  wavelengths  and  for  10-//m  (disks  only). 
Results  of  the  comparisons  are  tabulated  and  discussed. 

2.  Description  of  the  Problem 

As  a  part  of  the  CIAO  II  workshop,  members  of  the  development  community  were  invited  to 
participate  in  a  modeling  exercise  to  determine  the  agreement  among  the  various  models  in 
determining  basic  EM  properties  of  non-spherical  particles.  This  exercise  was  intended  to  serve 
as  a  baseline  study  from  which  to  progress  further.  Two  examples  were  chosen;  one  to  address 
IR  flakes  (disks)  and  one  for  finite  length  fibers.  All  fiber  calculations  were  performed  for  four 
wavelengths,  one  in  the  infrared  and  three  in  the  millimeter  wave  band.  Three  fiber  lengths  were 
chosen  to  span  the  dipole,  resonant,  and  diffraction  regimes.  In  all  cases,  the  fiber  was 
assumed  to  be  a  7-//m  diameter  cylinder  and  the  flakes  were  assumed  to  be  1 0-//m  diameter  by 
0.2-//m  thick  disks.  Calculations  for  the  flakes  were  performed  at  10-//m  wavelengths  only.  In 
all  cases,  the  calculations  were  made  assuming  an  ensemble  average  over  a  uniform  particle 
orientation  distribution  giving  what  is  referred  to  in  the  atmospheric  propagation  community  as 
a  tumble  average  of  all  quantities.  The  original  computational  plan  called  for  computing  all  four 
polarization  components  as  well  as  the  phase  function.  However,  for  model  output  comparison 
purposes,  we  focus  on  the  mass  extinction  coefficients  for  which  field  data  are  available  for 
comparison.  Details  of  the  specific  particle  characteristics  that  were  used  in  the  exercise  are 
summarized  in  table  C-1 .  Table  C-1  also  gives  values  for  the  bulk  conductivity  and  refractive 
indices,  which  happen  to  be  characteristics  of  graphite. 

Four  models  were  compared  in  the  millimeter  wave  regime  for  the  fiber  example;  AWAS, 
JP094M1 ,  IPHASE,  and  PA/V.  AWAS  is  based  upon  a  commercially  available  moment  method 
supplied  by  Bob  Richardson  of  the  JPO.  The  JP094M1  model  was  supplied  by  ARL  (Klett  and 
Sutherland)  and  is  a  modification  of  the  IPHASE  finite  cylinder  model  developed  by  Blair  Evans. 
The  IPHASE  model  is,  itself,  based  upon  the  variational  approach  published  over  the  years  in 
various  stages  of  development  by  Pederson,  Waterman,  and  Pederson  (PWP).  PA/V  is  also 
based  upon  the  PWP  approach  but  it  includes  a  special  current  function  that  eliminates  certain 
errors  found  in  the  earlier  versions  as  explained  by  Peter  Waterman  at  CIAO  I.  In  the  IR  regime 
(10-//m),  two  models  were  compared,  IPHASE  and  JP094S1,  which  are  both  based  upon  the 
infinite  cylinder  geometry  using  the  formal  solutions  published  in  the  literature  by  Bohren  and 
Huffman  (1983)  and  Kerker  (1969). 

Four  models  were  compared  for  flakes:  DDSCAT,  IPHASE1,  IPHASE2,  and  FLAKE.  The 
DDSCAT  model  is  based  on  the  Discrete  Dipole  Approximation  as  developed  by  Piotr  Flatau. 
The  IPHASE  models  were  supplied  by  Blair  Evans  and  represent  the  oblate  spheroid  and 
anomalous  diffraction  approximations,  respectively.  The  FLAKE  model  is  also  a  part  of  the 
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IPHASE  library.  At  the  workshop,  we  also  attempted  to  include  results  from  the  Discrete 
Green'?  Function  (DGF)  model  of  O’Brien  and  Goedecke  but  the  model  apparently  was  not  able 
to  handle  the  number  of  computations  required  for  the  particular  geometries  tested. 

3.  Results  for  Fibers  (Cylinders) 

For  the  10-A/m  wavelength,  IPHASE  and  JP094S1  computational  results  yielded  values  for  the 
mass  extinction  coefficient  (m^/g)  of  0.195  and  0.178,  respectively.  The  calculations  in  this 
region  should  be  reasonably  accurate  because  the  diffraction  regime  contains  numerous  valid 
approximation  methods  (including  the  infinite  cylinder). 

The  results  were  insensitive  to  the  actual  value  of  the  fiber  length  because  this  is  in  the  region 
of  applicability  for  the  infinite  cylinder  approximation.  The  difference  between  the  two  models, 
though  slight  (10  percent),  may  be  due  to  the  treatment  of  polarization  in  performing  the 
ensemble  averaging.  This  same  problem  arose  in  discussions  involving  the  phase  function. 
Several  ideas  presented  at  the  workshop  are  currently  being  investigated  to  resolve  the 
differences. 

Results  for  the  millimeter  wave  regime  are  summarized  in  table  C-2  where  all  four  applicable 
models  were  used.  The  mass  extinction  coefficient  is  the  quantity  compared  at  three  millimeter 
wavelengths.  For  all  table  C-2  cases  the  mass  extinction  coefficient  was  determined  from  the 
calculated  extinction  cross-section  and  a  bulk  mass  density  of  1 .78  g/cm^  which  is  appropriate 
for  graphite. 

From  inspection  of  the  first  three  columns  of  table  C-2,  which  correspond  to  the  shortest  fiber 
length  (3.175mm),  the  results  for  the  three  variational  models  (JP094M1,  IPHASE2,  PNSl,  and 
AWAS)  are  within  10  percent  agreement.  Earlier  results  for  the  IPHASE  calculations  as 
presented  at  the  workshop  were  later  found  to  be  in  error  by  a  factor  of  /r/4  sin(0)  and  are  shown 
as  a  separate  model  that  is  identified  as  IPHASE1.  The  error  was  caused  by  an  oversight  in 
defining  the  solid  angle  for  integration  over  the  particle  orientation.  The  results  for  the  IPHASE1 
model  are  listed  for  the  sake  of  completeness  because  they  were  presented  at  the  CIAO  II 
workshop  but  they  should  othenwise  be  ignored.  Previously  reported  values  for  AWAS  were  also 
later  found  to  be  in  error  and  have  been  corrected  here. 

The  results  for  the  intermediate  fiber  length  (6.35-mm)  show  a  more  marked  disparity  among  the 
three  variational  models  with  the  P/W  results  being  significantly  higher  than  the  JP094M1  or 
IPHASE2,  at  the  shortest  wavelength.  Also  the  AWAS  results,  are  generally  in  good  agreement 
with  all  others.  In  cases  where  there  is  marked  disparity  among  methods,  the  JP094M1  and 
IPHASE2  results  tend  to  agree  with  each  other  while  the  AWAS  Method  of  Moments  based  code 
calculations  tend  to  agree  more  closely  with  the  Pedersen-Waterman  results  which  use  an 
improved  current  function.  It  is  worth  noting  that,  overall,  the  results  for  the  two  shorter 
wavelengths  are  in  fair  agreement  with  field  data  that  estimate  values  of  around  2  to  2.5  m%. 

Results  for  the  longest  fiber  length  (12.75-mm)  show  even  greater  disparities  between  the 
variational  methods  with  the  P/W  results  being  over  a  factor  of  two  higher  than  the  JP094M1  or 
IPHASE2.  This  is  due  to  the  improved  treatment  of  the  fiber  current  function  employed  by 
Waterman  in  the  P/Wformuation.  The  original  version  presented  errors  in  certain  wavelength 
regimes.  However,  the  cause  of  this  disparity  and  that  with  AWAS  is  still  being  investigated. 
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A  word  about  the  material  parameter  model  for  this  calculation  exercise  is  in  order.  For  all  of 
the  microwave  and  mm-waye  calculations,  the  carbon  fiber  electrical  parameters  were  assumed 
to  follow  a  Debye  dielectric  model  rather  than  the  Drude  model  where  conductivity  is  essentially 
determined  by  its  DC  value  and  dielectric  displacement  current  is  negligible.  While  this  is 
appropriate  for  comparison  of  different  calculation  tools,  there  is  some  doubt  that  the  Debye 
model  is  the  most  appropriate  for  describing  carbon  fiber  electrical  behavior  in  the  microwave 
and  mm-wave  regime.  Differences  between  the  two  models  become  more  severe  as  one  goes 
lower  in  frequency.  In  other  calculation  and  measurement  exercises,  the  substitution  of  DC 
conductivity  rather  that  Debye  model  based  parameters  appeared  to  yield  results  which  were 
more  nearly  in  agreement  with  measurement. 

As  a  matter  of  completeness,  table  C-3  shows  the  results  of  some  intermediate  calculations 
supplied  by  Blair  Evans  using  the  various  modules  in  the  IPHASE  library.  Table  C-3  includes 
results  for  the  infinite  cylinder,  the  finite  cylinder,  and  the  prolate  spheroid  covering  all  the  cases 
suggested  in  the  exercise  where  the  various  approaches  are  assumed  valid. 

4.  Results  for  Flakes  (Disks) 

Results  for  IR  (10-A/m)  flakes  giving  values  for  the  mass  extinction,  absorption,  and  scattering 
coefficients  (m^/g)  using  the  DDSCAT,  IPHASE,  and  FLAKE  approximations  are  summarized  in 
table  C-4. 

Note  that  in  table  C-4  the  DDSCAT  method  yielded  a  value  of  2.40  for  the  scattering  component 
and  1 .08  for  the  absorption  component  giving  an  albedo  of  0.43  as  compared  to  a  value  of  0.3 
for  the  FLAKE  model.  Currently,  it  is  not  clear  why  all  methods  are  net  in  closer  agreement. 
It  is  interesting  to  note  that  field  measurements  yield  values  ranging  between  1  and  3  m^/g, 
which  are  not  signficantly  greater  than  the  spread  in  values  in  table  C-4  if  FLAKE  results  are 
excluded.  An  additional  factor  affecting  the  comparison  of  theoretical  and  experimental  data  are 
the  validity  of  assumptions  of  disk  size  and  shape  for  the  IR  flake  materials. 

Intermediate  data  showing  plots  of  the  extinction,  absorption,  and  scattering  cross-sections  and 
phase  function  asymmetry  factor  as  a  function  of  orientation  angle  are  shown  in  figure  C-1. 
These  plots  were  supplied  by  Pitor  Flatau  using  data  based  upon  the  DDSCAT  model,  which  is 
believed  to  be  the  most  accurate  representation  of  any  of  the  models  presented  here.  Even  so, 
the  DDSCAT  mass  extinction  coefficient  is  about  20  percent  higher  than  field  measurements  of 
graphite  flakes. 

The  phase  function  plots  of  figure  C-1  appear  well  behaved,  beginning  at  some  initial  high  value 
near  zero  (broadside)  and  diminish  as  the  angle  increases.  The  assymetry  parameter  varies 
from  near  zero  (isotropic  scattering)  at  the  lower  angles  to  near  unity  (strong  fonward  scattering) 
at  the  higher  angles.  Plots  of  the  modeled  phase  functions  for  the  four  polarization  cases  are 
plotted  in  figure  C-2. 
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Table  C-1.  Input  data  for  the  modeling  exercise 


Incident 

Radiation  Particle  Characteristics 

Circular 

Disk  Fibers  Bulk  Admittance  mho/m 

Complex 

Index  of  Refraction 

Wavelength  Diameter;  lOjitm 

Thickness:  0.2/i  m  Diameter;  7/im 

Lengths; 

3.175,  6.35,  and  12.7mm  a  jwe 

n 

k 

lOlim  X  X  3.47  X  10^ 

4.791 

3.513 

3.19mm  NA  X  7.67x10“ 

2.45x10“ 

101 

72.77 

8.57mm  NA  X  4.87x10“ 

4.18x10“ 

165.7 

75.69 

18.75mm  NA  X  1.87x10“ 

3.58x10“ 

205.8 

51.2 

Table  C-2.  Results  of  the  CIAO  II  workshop  modeling  exercise  giving  mass  extinction  coefficients 
for  fiber  lei^hs  of  3.175-mm,  6.35-mm,  and  12.75-mm. 


- PREDICTED  MASS  EXTINCTION  COEFFICIENTS,  m"/g 

MODEL 

1  = 

=  3.175nim 

1  = 

6.35mm 

1  = 

12.75mm 

X, 

^2 

X3 

X, 

^2 

X3 

X, 

X2 

^3 

JP094M1 

2.56 

1.74 

0.206 

1.65 

3.63 

0.80 

1.36 

1.74 

1.31 

IPHASE2 

2.54 

1.73 

0.222 

1.65 

3.60 

0.80 

1.35 

1.73 

1.30 

P/W 

2.72 

1.73 

0.245 

2.51 

3.58 

0.81 

2.45 

4.41 

1.31 

AWAS2 

2.32 

1.56 

0.227 

2.16 

3.22 

0.76 

2.09 

3.85 

1.25 

AWASl 

1.63 

3.02 

2.23 

IPHASEl 

3.90 

2.57 

0.340 

2.66 

5.30 

1.21 

2.20 

2.82 

1.94 

X,  =  3.19mm  (94GHz);  Xj  =  8.57mm  (35GHz);  X3  =  18.75mm  (16Ghz) 
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Table  C3.  Intermediate  results  supplied  by  Blair  Evans;  Lj  =  3.175mm,  Lj  =  6.350mm, 
L3  =  12.75mm 


Model 

Qex. 

Qsca 

Qabs 

a 

10 

^im 

3.19 

mm 

8.57 

mm 

18.75 

mm 

10 

^m 

3.19 

mm 

8.57 

mm 

18.75 

mm 

10 

^m 

3.19 

mm 

8.57 

mm 

18.57 

mm 

10 

|im 

3.19 

mm 

8.57 

mm 

18.75 

mm 

H 

2.6977 

58.274 

112.81 

43.180 

1.9919 

15.818 

16.357 

3.1897 

.70578 

42.456 

96.451 

39.991 

.22276 

4.8120 

9.315 

3.5655 

2.0326 

.00238 

1.2792 

1,96'^ 

j^^jj 

9.46''^ 

.75341 

.00238 

■ 

1.63'^ 

.1678 

Finite 

Cylindei 

L, 

.41176 

30.757 

20.929 

2.688 

.35068 

7.345 

1.418 

11 

.00610 

23.411 

19.511 

2.688 

.02944 

2.543 

1.73 

.2222 

hi 

.49937 

19.096 

43.572 

9.782 

.48865 

4.5137 

5.064 

.2432 

.01082 

15.396 

38.508 

9.539 

.03570 

1.645 

3.600 

0.8082 

L3 

.49221 

16.377 

20.921 

15.769 

.48221 

3.888 

2.124 

0.8456 

.00925 

12.490 

18.797 

14.923 

.03511 

1.353 

1.328 

1.307 

Prolate 

Spheroid 

L, 

2.4698 

2.1950 

1.2158 

0.4651 

■ 

■ 

■ 

- 

■ 

■ 

■ 

.24154 

.21354 

.11945 

.04519 

hi 

2.4700 

2.1950 

1.2158 

0.4651 

- 

- 

- 

- 

- 

- 

* 

.24154 

.21354 

.11945 

.04519 

L3 

2.4700 

2.1950 

1.2158 

0.4651 

> 

- 

- 

- 

- 

- 

- 

• 

.24154 

.21354 

.11945 

Table  C-4.  Summary  of  results  for  the  disk  modeling  exercise 


Model 

aabs(mVg) 

asca(mVg) 

OtexXmVg) 

DDSCAT 

1.08 

2.40 

3.48 

IPHASE  1 

2.30 

IPHASE  2 

2.28 

FLAKE 

9.22 

3.88 

13.1 
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Figure  C-1.  Intermediate  data  for  the  flake  example  showing  cross-sections  and 
asymmetry  parameter  as  a  function  disk  orientation. 
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Appendix  D 


Review  of  Existing 
Measurement  Techniques 


76 


1.  Purpose  and  Scope 

This  appendix  describes  experimental  data  and  data  analyses  methods  routinely  used  for 
evaluating  obscuration  models.  Most  of  the  data  rom  field  experiments  that  can  be  used  to 
evaluate  the  accuracy  of  electromagnetic  wave  propagation  model  prediction  are  multispectral 
extinction.  Effective  electromagnetic  wave  propagation  models  must  be  able  to  predict  this  basic 
integral  parameter.  Multispectral  extinction  is  used  to  estimate  obscurant  logistics  requirements 
and  to  define  the  operational  threshold  of  obscured  sensors.  Therefore,  it  is  important  to  have 
a  clear  understanding  of  how  multispectral  extinction  is  normally  analyzed  and  used  to  test 
propagation  models. 

The  largest  data  base  of  field  measurements  available  for  model  prediction  comparisons  is  that 
for  multispectral  extinction.  Experimental  techniques  for  acquiring  these  data  are  reviewed,  and 
methods  of  data  analyses  are  discussed.  It  is  shown  how  attenuation  data  are  converted  to  mass 
extinction  coefficients.  Sources  of  experimental  data  are  described  and  laboratories  and 
agencies  that  have  contributed  significantly  to  obscurant  data  are  discussed.  The  application 
of  the  experimental  data  base  and  measurement  capability  to  validation  of  predictive  models  is 
described  for  the  scattered  electromagnetic  vector  field  as  a  function  of  particle  orientation,  the 
orientation  distribution  of  fibers  in  the  atmosphere,  particle  number  density  after  dissemination, 
the  dielectric  constant  of  obscurant  materials,  the  obscurant  characteristic  polarization,  and  error 
models  of  instrumentation  used  to  gather  the  data. 

2.  Multispectral  Extinction  Fieid  Measurements 

Figure  D-1  shows  a  typical  experimental  arrangement  for  the  measurement  of  multispectral 
extinction  through  an  atmospheric  obscurant.  The  obscurant  aerosol  is  released  upwind  from  a 
multispectral  transmissometer  line  of  sight.  Typically,  transmittance  through  the  aerosol  as  it 
drifts  through  the  line  of  sight  is  sampled  at  1  to  10  Hz  rates.  The  data  are  acquired  as  a  function 
of  time  from  some  arbitrary  zero  reference.  Each  data  point  represents  the  average  of  a  number 
of  sampled  points  taken  at  a  much  higher  rate  than  the  sampling  frequency. 

The  tranmissometer  spectral  range  may  be  covered  by  a  continuous  broad  band  or  by  a  series 
of  narrow  bands.  Narrow  band  measurements  are  needed  when  the  spectral  extinction  varies 
significantly  within  the  measured  band.  Broad  band  measurements  are  sufficient  when  the 
spectral  extinction  is  constant  within  the  measured  band. 

Path  lengths  through  obscurants  are  typically  10  to  100  m  long.  Line  of  sight  path  lengths 
(distance  between  transmitter  and  receiver)  range  between  500  and  3000  m.  For  visual  to 
infrared  bands,  field  test  systems  measure  at  0.4  to  0.7, 1 .06,  1 .54,  3  to  5,  8  to  14,  and  10.6//m 
wavelengths.  At  millimeter  wavelengths,  transmittances  are  obtained  at  35  and  95  GHz. 

To  determine  the  mass  extinction  coefficient,  the  line  of  sight  concentration  length  must  be 
measured  simultaneously  with  the  transmittances.  This  is  accomplished  with  an  array  of  point 
mass  concentration  sensors  or  by  transmittance  at  a  wavelength  for  which  the  mass  extinction 
coefficient  is  known.  Point  samplers  are  commonly  used  to  measure  visual  to  infrared  obscurant 
mass  concentrations.  These  devices  are  also  used  for  measuring  number  density  of  millimeter 
wave  obscurants. 
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Figure  D-1.  Typical  transmissometer  system  in  a  field  test  scenario  for  the  measurement  of 
multispectral  extinction. 
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Errors  in  the  concentration  length  measured  with  a  sensor  array  arise  because  of 
inhomogeneities  in  the  obscurant  concentration  distribution  between  samplers.  These 
concentration  length  errors  vary  in  ‘nagnitude  from  a  few  percent  to  factors  of  two. 
Concentration  lengths  measured  with  transmissometers  are  mac;  for  millimeter  wave 
obscurants.  In  the  case  of  millimeter  wave  obscurants,  the  mass  extinction  coefficient  can  be 
estimated  within  a  few  percent  relative  uncertainty  using  the  geometric  optical  inversion  (GOI) 
technique.  The  GOI  technique  works  in  this  case  because  the  particles  are  of  a  size  that  the 
extinction  efficiency  is  approximately  2  at  visual  to  infrared  wavelengths. 

3.  Multispectral  Extinction  Data  Anaiyses 

For  a  well  collimated  transmissometer  where  data  are  produced  that  have  no  multi-scatter 
artifacts,  the  Beer-Bouguer  transmission  law  can  be  applied  to  the  data  analyses.  Minor 
manipulation  of  equation  (1)  in  section  3,  Volume  I,  shows  that  the  Beer-Bouguer  law 
transmission  description  for  a  homogeneously  distributed  obscurant  material  is 

r=exp(-opy/) 

D-1 


where 

T  =  the  ratio  of  transmitted  to  incident  irradiance 

o  =  the  mean  extinction  cross  section  of  the  attenuating  material 
=  the  aerosol  particle  number  density 
L  =  the  path  length 

In  equation,  D-1,  the  extinction  cross-section  is  a  function  of  wavelength.  The  wavelength 
dependence  is  assumed  in  the  following  equations  and  not  written  explicitly. 

Aerosol  particle  number  densities  are  difficult  to  measure  experimentally:  therefore,  it  is  common 
practice  to  transform  the  variables  into  those  related  to  mass  (which  is  straightforward  to 
measure).  The  mass  extinction  coefficient  a,  is  the  mean  extinction  cross-section  per  mean 
particle  mass.  For  spherical  particles  the  mass  extinction  coefficient  is 


TZ 

gPoF3 


(D-2) 


where 

Po  =  the  material  density  of  the  obscurant 

P3=  the  third  moment  of  the  particle  size  distribution 
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The  spherical  particle  mass  concentration  is 


(D-3) 

Using  equations  (D-2)  and  (D-3)  in  equation  (D-1)  the  expression  normally  used  in  obscurant 
analyses  is  as  follows: 


7=exp(-aCZ.) 

(D-4) 

The  concentration  length  (CL),  as  shown  in  equation  (D-1),  is,  in  fact,  an  integral  expression 
because  the  concentration  is  not  homogeneously  distributed  along  the  transmittance  line  of  sight. 
CL  is 


L 

CL=fC{r)dl 

0 

(D-5) 

Field  measurements  produce  values  of  T  and  CL.  The  mass  extinction  coefficient,  predicted 
from  electromagnetic  wave  propagation  models,  is  from  experimental  data: 


(D-6) 

Although  equations  (D-1)  to  (D-6)  are  relatively  simple,  there  are  a  number  of  subtle  points  that 
must  be  considered  to  effectively  compare  experimental  results  with  theoretical  predictions.  First 
consider  the  impact  of  error  in  transmittance  and  concentration  length  on  the  relative  uncertainty 
in  mass  extinction  coefficient.  The  root-mean-square  (RMS)  relative  uncertainty  in  mass 
extinction  coefficient  is 


a  ; 


■-1 


dVTf(HCL)\^ 
ln(7)j  [{CL) 


(D-7) 


where  6x/x  (x  =  a)  is  the  relative  uncertainty  in  the  variable,  and  a  bar  signifies  a  mean  value. 
Figure  D-2  plots  the  relative  error  in  mass  extinction  coefficient  as  a  function  of  transmittance 
and  selected  relative  uncertainties  in  transmittance  and  CL. 


Figure  D-2  shows  that  the  relative  uncertainty  varies  with  the  transmittance  and  has  a  minimum 
value  when  the  transmittance  is  about  0.37.  The  relative  uncertainty  in  mass  extinction 
coefficient  increases  rapidly  near  transmittances  of  1  and  when  the  transmittance  is  less  than 
about  0.1.  Therefore,  the  relative  uncertainty  in  mass  extinction  must  be  determined  by 
integrating  equation  (D-7)  over  the  distribution  of  transmittance  values  obtained  for  the  range  of 
CL  used  in  the  estimate. 
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TRANSMITTANCE 


Figure  D-2.  Relative  error  in  measured  Mass  extinction  coefficient  for  05  uncertainty  in 
transmittance  and  selected  relative  uncertainties  in  concentration  length. 
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The  second  point  to  consider  in  evaluating  the  mass  extinction  coefficient  is  its  dependence  on 
the  obscurant  particle  size  distribution.  The  shape  of  the  particle  size  distribution  defines  the 
third  moment  //j  of  the  size  distribution  and  is  used  to  compute  the  mean  value  of  the  extinction 
cross-section  in  equation  D-2.  It  is  possible  to  experimentally  determine  the  mass  extinction 
coefficient  through  a  measurement  of  the  particle  size  distribution  if  the  particle  material  density 
is  known  or  assumed.  However,  particle  size  distribution  measurements  sample  a  very  small 
volume  of  air  relative  to  that  of  a  transmissometer,  and  the  size  distribution  measurement  itself 
is  very  sensitive  to  sampling  and  measurement  errors. 

The  analyses  given  by  equations  (D-1)  to  (D-7)  is  necessary  if  absolute  values  of  mass  extinction 
coefficient  are  to  be  determined.  The  uncertainties  in  CL  measurements  of  many  obscurants  or 
particle  size  distribution  measurements  make  the  determination  of  absolute  values  of  mass 
extinction  coefficient  difficult  and  these  determinations  often  have  a  high  level  of  uncertainty. 
Thus,  methods  using  linear  regressions  of  multispectral  attenuation  measurements  have  been 
developed  that  eliminate  the  need  to  directly  measure  concentration  lengths.  The  following 
discussion  outlines  this  method  of  determining  relative  mass  extinction  coefficients. 

Assume  that  a  transmissometer  is  capable  of  simultaneously  measuring  two  or  more  spectral 
transmittances.  The  negative  natural  logarithm  of  a  spectral  transmittance  is  called  the  optical 
depth.  The  ratio  of  optical  depths  for  two  simultaneously  measured  transmittances  is 

-ln(7i:Xi))  tt{X^)CL 
aiX^CL 

(D-8) 

Assume  that  the  optical  depth  for  is  the  dependent  variable.  Equation  (D-8)  can  then  be 
written  as  a  linear  equation  of  the  form 


y=mx+b 


(D-9) 


where 


y=-ln(7lX0) 

x=-\nmx^) 

«(^i) 

a{X^ 

b=0 

Much  information  is  gained  by  plotting  y  versus  x  and  performing  a  linear  regression.  First,  if 
there  is,  in  fact,  in  a  linear  relationship,  then  the  Beer-Bouguer  transmission  law  applies  to  the 
data.  Second,  the  correlation  coefficient  from  the  regression  is  a  measure  of  the  noise  and 
scatter  of  the  data  between  the  two  spectral  bands.  The  correlation  coefficient  provides  a 
measure  of  data  quality  between  the  measurements.  Third,  if  b  0,  then  there  is  clear 
indication  that  the  transmissometer  is  operating  incorrectly  or  that  data  artifacts  should  be 
expected.  Fourth,  the  slope  from  the  regression  is  the  ratio  of  mass  extinction  coefficients  for 
the  two  spectral  bands.  Figure  D-3  illustrates  a  typical  regression  plot  of  data  measured  in  a 
field  experiment.  Information  derivable  from  these  type  of  analyses  is  of  sufficient  value  that 
virtually  all  field  measured  transmittance  data  are  now  reviewed  using  this  technique. 
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Evaluation  of  particle  mass  concentrations  or  number  densities  for  millimeter  wave  and 
microwave  obscurants  in  field  tests  is  accomplished  by  measurements  at  an  array  of  points  along 
a  line  of  sight  using  optical  nephelometers  calibrated  to  count  single  particles  or  to  measure 
large  numbers  simultaneously  for  mass  concentrations.  Mean  mass  or  number  concentrations 
along  a  propagation  path  can  be  computed  using  the  Geometric  Optical  Inversion  (GOI) 
technique  on  optical  (visible  or  infrared)  transmittance  data.  The  GOI  can  be  used  for  millimeter 
wave  or  micrswave  obscurants  when  it  is  assumed  the  particles  are  much  larger  than  an  optical 
wavelength  and  are  randomly  oriented  (see  Wright  and  Butters  in  the  Bibliography  in 
appendix  B).  Under  these  assumptions,  the  mass  extinction  coefficient  (extinction  cross-section 
per  particle  mass)  can  be  accurately  estimated.  It  is 


a  = 


2 

Po^ 


(D-10) 


where  is  the  particle  material  density  and  D  is  the  diameter.  Equation  (D-10)  shows  that  the 
mass  extinction  coefficient  is  independent  of  incident  wavelength  and  particle  length  when  D»/t. 
Equation  (D-1 0)  can  be  used  in  the  Beer-Bouguer  transmission  law  to  obtain  CL.  The  result  is 


CL=  ^ 

(D-11) 

The  weaknesses  in  estimating  CL  represented  by  equation  (D-11)  are  1)  the  particles  must  be 
randomly  oriented,  and  2)  the  transmission  line  of  sight  must  be  representative  of  the  line  of 
sight  of  the  millimeter  or  microwave  path.  Furthermore,  if  number  density  is  to  be  inferred  from 
the  CL  then  absolute  path  length  through  the  obscurant  is  required  and  particle  length  must  be 
known.  Applications  of  the  GOI  technique  to  chamber  measurements  of  millimeter  wave  and 
microwave  obscurants  are  straightforward.  Path  lengths  in  the  chamber  are  known.  Turbulence 
in  the  chamber  with  mixing  fans  should  assure  random  orientation.  Particle  releases  in  the 
chamber  are  sufficiently  gentle  so  that  particles  are  not  broken  and  the  size  distributions  are 
monodisperse  in  length  and  the  particles  are  not  agglomerated. 

Uncertainties  in  the  GOI  technique  include  those  determined  by  non-random  particle  orientation, 
particle  length  distributions,  and  particle  agglomerations.  Data  for  particle  size  distributions  of 
millimeter  wave  and  microwave  obscurants  are  extremely  limited.  Analyses  of  transmittance  data 
for  field  dispersions  of  these  obscurants  strongly  suggest  variations  in  size  distributions  result 
from  variations  in  length  (resulting  from  broken  particles)  and  variations  in  diameter 
(agglomerated  particles).  Data  reflecting  direct  measurement  of  these  parameters  in  field 
measurements  is  sparse  or,  in  some  cases,  nonexistent. 
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Figure  D-3.  Example  regression  plot  of  multispectral  optical  depths  to  determine  relative 
mass  extinction  coefficients. 
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4.  Radar  Measurements 


If  the  transmissometer  is  replaced  by  a  radar  such  that  scattering  vs.  range  may  be  obtained, 
in  addition  to  transmission,  more  complete  information  about  the  mass  concentration  in  the  cloud 
and  a  more  complete  overall  picture  to  the  dissemination  may  be  obtained.  However,  radar 
measurement  is  more  complicated  than  transmissometry.  Issues  involved  in  determining 
whether  simple  transmissometry  or  more  complicated  scattering  measurement  is  required  along 
with  an  overview  of  scattering  measurement  is  described  below. 

Typically,  the  radar  should  be  capable  of  range  resolution  on  the  order  of  0.5  to  2  meter  in  order 
to  determine  the  mass  distribution  in  the  cloud  with  suitable  definition  and  also  to  aid  in 
correcting  for  absorption  of  the  signal  with  the  cloud  to  determine  the  true  scattering.  Workers 
form  NSWCDD  (Richardson,  1995,  Bean,  1993)  have  shown  that  these  requirements  may  be 
easily  met  with  FMCW  radar  architecture.  Workers  from  GTRI  and  Eglin  AFB  (Perry,  1992, 
Mijangos,  1991)  have  performed  similar  measurements  using  a  pulsed  radar  with  a  logarithmic 
receive.  Figure  D-4  shows  an  elevation  view  of  a  typical  field-test  installation.  Several  corner 
reflectors  are  placed  down  range  so  that  the  transmission  and  scattering  are  obtained 
simultaneously.  The  scattering  measurement  provides  a  2-d  map  of  material  concentration. 
Spatial  resolution  is  determined  by  the  range  resolution  along  the  beam  and  the  beamwidth 
along  the  direction  of  cloud  motion.  By  using  several  corner  reflectors  at  different  elevations,  an 
estimate  (neglecting  diffraction  effects)  of  the  vertical  extent  of  the  cloud  may  also  be  obtained. 
A  typical  A-scope  display  is  shown  in  Figure  D-5.  When  the  cloud  passes  through  the  beam,  its 
location  along  the  beam  and  its  extent  are  readily  visible  from  the  scattering  return  while  the 
attenuation  is  observed  directly  by  the  reduction  in  signal  level  from  the  comer  reflector  peaks, 
the  beamwidth  of  the  radar  beam  must  be  carefully  characterized  so  that  the  spot  size  is  a 
known  function  of  range.  Further,  it  is  desirable  that  the  beam  be  narrow  enough  so  that  ground 
reflection  effects  and  ground  clutter  are  negligible.  Ideally,  the  center  of  the  beam  should  be  at 
the  same  elevation  as  the  nephelometer  line. 

Interestingly,  the  data  obtainable  from  scattering  measurement  is  complementary  to 
nephelometer  data.  If  one  requires  good  spatial  resolution  of  the  material  distribution  in  a  cloud, 
and  relies  solely  on  nephelometers,  a  great  many  nephelometers  are  required  since 
nephelometer  lines  are  also  required  to  be  quite  long.  On  the  other  hand,  if  one  relies  on  the 
range  resolution  scattering  for  cloud  material  distribution  measurement,  the  nephelometers  may 
be  relied  upon  for  more  normalization  data  to  check  for  consistency  of  results.  High  range 
resolution  and  good  normalization  may  be  obtained  without  an  infinite  number  of  nephelometers. 

With  a  suitably  calibrated  radar,  (and  corrections  for  attenuation  within  the  cloud)  one  measures 
RCS  per  unit  volume  of  cloud  material.  The  pulse  volume  is  determined  by  the  spot  size  and 
the  range  resolution  which  is  determined  by  the  FFT  block  size  and  windowing  function  chosen 
for  processing  the  FMCW  radar  burst. 

Because  of  clumping  effects,  scattering  measurement  provides  a  better  estimate  of  mass 
concentration  in  a  cloud  and  a  much  better  overall  picture  of  the  dissemination  than  absorption 
measurement  alone,  particularly  for  resistive  obscurant  materials.  It  is  well  known  that  materials 
do  not  always  disseminate  as  single  fibers  but  often  clump  together  into  multiple-fiber  bundles 
where  several  (n)  fibers  are  positioned  together  in  a  cord  wood  fashion.  For  typical  parameter 
values  clumping  has  profound  influence  on  the  absorption  coefficient  as  shown  in  Figure  D-6  but 
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Figure  D-4.  Elevation  (side)  view  sketch  of  a  typical  radar  field  set-up, 
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Figure  D-6.  Calculated  a  and  a,  (mVgm)  verses  number  of  fibers  in  a  bundle. 
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not  so  much  influence  on  the  scattering  coefficient.  In  a  typical  dissemination  event  there  is  a 
"spectrum"  of  chaff  bundles  with  various  "n"  values.  An  important  parameter  in  determining 
absorption  coefficient,  is  the  resistance  per  unit  length  of  the  fiber  bundle.  Assuming  that  the 
material  is  optimized  for  single  fiber  dissemination,  a  multifiber  bundle  weighs  more  and  has  a 
lower  than  optimal  resistance  per  unit  length  to  provide  attenuation. 

In  a  practical  field  measurement  or  even  in  a  chamber,  one  does  not  have  a  very  good  way  of 
determining  the  average  number  of  fibers  per  bundle.  Counting  the  number  of  fibers  in  a  bundle 
under  a  microscope  is  very  laborious.  The  calculated  results  of  Figures  D-6  show  that  the 
scattering  cross  section  per  unit  mass  has  a  much  weaker  dependence  on  the  number  of  fiber, 
i.e.  the  clumping  number  (n)  than  does  the  absorption.  Since  it  is  laborious  to  get  a  good 
estimate  of  n  from  microscopy,  mass  concentration  should  be  based  on  corrected  sea  taring  data 
rather  than  attenuation  data. 

With  a  suitably  calibrated  radar,  one  measures  volume  scattering  coefficient  I’m^/m^)  and  then 
corrects  the  measured  scattering  for  attenuation  such  that  the  attenuation  at  ;  !  ;8  far  side  of  the 
cloud  matches  that  attenuation  observed  on  the  corner  reflectors.  In  performing  this  correction, 
one  naturally  arrives  at  the  ratio  of  scattering  and  extinction  cross  sections.  The  data  of  Figure 
D-6  shows  that  this  ratio  is  clearly  a  function  of  how  much  clumping  is  present.  Reflectivity  or 
"blackness"  of  a  cloud  which  depends  on  the  ratio  of  scattering  and  absorbing  cross  sections  is 
clearly  dependent  on  the  amount  of  clumping  in  the  cloud. 

In  material  development  studies  where  the  reflectivity  of  the  material  is  of  interest  or  in 
measurements  to  evaluate  the  performance  of  cloud  formations  models  where  mass  per  unit 
volume  is  of  interest  it  is  clearly  important  to  perform  scattering  measurements  in  addition  to  the 
more  usual  attenuation  measurements. 

Uncertainties  in  the  GOI  technique  include  those  determined  by  non-random  particle  orientation, 
particle  length  distributions,  and  particle  agglomerations.  Data  for  particle  size  distributions  of 
millimeter  wave  and  microwave  obscurants  are  extremely  limited.  Analysis  of  transmittance  data 
for  field  dispersions  of  these  obscurants  strongly  suggests  variations  in  size  distributions  that 
result  from  variations  in  length  (resulting  from  broken  particles)  and  variations  in  diameter 
(agglomerated  particles).  Data  reflecting  direct  measurement  of  these  parameters  in  field 
measurements  are  sparse  or  in  some  cases  nonexistent. 

5.  Data  Sources 

Sources  of  experimental  data  to  verify  electromagnetic  wave  propagation  models  fall  into  two 
broad  categories:  1)  laboratory  measurements,  and  2)  field  trials.  Data  from  both  sources 
can  be  grouped  into  three  categories:  1)  material  characteristics,  2)  electromagnetic  wave 
response  characteristics,  and  3)  figures  of  merit.  Figures  of  merit  are  derived  from  material 
characteristics  and  electromagnetic  wave  responses.  Table  D-1  provides  a  summary  of  the 
sources  of  data  in  these  categories.  The  agencies  listed  in  table  D-1  either  acquired  the  data 
using  government  facilities  or  through  contractors  sponsored  by  the  agency. 
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Table  D-1.  Summary  of  obscurant  data  sources  for  evaluation  of  electromagnetic  wave 
propagation  models. 
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[. . . . 
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CANADA 
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MATERIAL  CHARACTERISITICS 

m 

ELECTRICAL  CONDUCTIVTIY 

DELECTRIC  CONSTANTS 

B 

B 

B 

REFRACTIVE  INDICES 

B 

CHEMICAL  COMPOSITION 

B 

5 

CRYSTALUNE  STRUCTURE 

6 

DENSITY 

X 

7 

PARTICLE  SHAPE 

B 

B 

PARTICLE  SIZE,  SIZE  DISTRBtrnON 

X 

X 

X 

X 

X 

X 

B 

PARTICLE  NUMBER  DENSITY 

ELECTROMAGNETIC  WAVE  RESPONSE 

S11  PHASE  FUNCTION 

X 

X 

■iS 

MUELLER  MATRIX  PHASE  FUNCTIONS 

BIS 

X 

B 

EXTINCTION  CROSS-SECTIONS 

HH 

B 

X 

B 

X 

X 

RADAR  CROSS-SECTION 

mm 

B 

X 

X 

FIGURES  OF  MERIT 

VOLUME  EXTINCTION  COEFFICIENT 

B 

mi 

B 

X 

BiS 

MASS  EXTINCTION  COEFFIOENT 

B 

B 

B 

B 

X 

MM 

MASS  RADAR  COEFFICIENT 

X 

B 

X 

DRAG  COEFFICIENTS 

PARTICLE  ORENTATION 

B 

B 

X 

X 

CROSS-SECTION  RATIO 

B 

IBI 

FOURER  POWER  SPECTRUM 

B 

*Note:  HARC  has  not  been  used  for  obscurants. 
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5.1  Facilities 


Data  acquired  by  the  U.S.  Army  Research  Laboratory  (ARL),  Battlefield  Environment  Directorate 
(BED)  are,  for  the  most  part  multispectral  transmittance,  stereo  imagery  for  determining 
obscurant  transport  and  diffusion,  point  measurement  of  mass  concentration,  particie  size 
distribution,  and  meteorological  measurements  in  support  of  Smoke  Week  test  programs  and 
other  electro-optical  sensor  evaluation  programs. 

Dugway  Proving  Grounds  (DPG)  has  obtained  multispectral  transmittance,  point  concentration 
and  dosage  measurements  for  test  programs  at  DPG  and  for  other  electro-opticai  sensor  test 
programs,  such  as  the  Smoke  Weeks. 

Edgewood  Research,  Development,  and  Engineering  Center  (ERDEC)  has  generated  voluminous 
obscurant  data  through  laboratory  chamber  measurements,  small  scale  atmospheric  tests,  and 
major  test  programs  such  as  Smoke  Weeks.  ERDEC  has  sponsored  research  for  measurements 
of  material  characteristics  and  has  developed  a  small  nephelometer  capable  of  Mueller  matrix 
measurements  at  visible  laser  wavelengths. 

The  Air  Force  represented  by  Wright-Patterson  Air  Force  Base  (AFB)  and  Eglin  AFB  has 
conducted  numerous  chaff  (scattering  cylinders)  tests  in  anechoic  chambers  and  under  flight 
conditions.  Eglin  AFB  has  provided  millimeter  wave  radar  measurements  of  absorbing  cylinder 
tests  conducted  during  Smoke  Week  tests. 

The  Institute  for  Space  Science  and  Technology  (ISST),  Space  Astronomy  Laboratory, 
Gainesville,  FL,  operates  a  3.1 -cm  wavelength  microwave  nephelometer  capable  of  measuring 
phase  functions  for  irregularly  shaped  particles  scaled  for  extinction  efficiency  measurements  in 
the  resonant  region. 

ISST  also  measures  the  dielectric  constants  of  the  materials  used  to  make  the  particle  models 
used  in  the  measurements.  Phase  functions  measured  in  this  laboratory  have  provided  detailed 
data  for  evaluation  of  theoretical  models  of  resonant  region  scattering. 

The  Houston  Area  Research  Consortium  (HARC)  have  obtained  data  similar  to  that  of  ISST. 
However,  the  HARC  emphasis  has  been  on  measuring  bistatic  radar  cross-sections  of  object 
shapes  of  interest  to  the  military  and  has  not  measured  obscurants. 

Data  acquired  in  the  United  Kingdom  include  both  laboratory  and  field  measurements  with 
multispectral  transmissometers  measuring  visible  and  emissive  infrared  obscurants,  and  radar 
systems  observing  chaff  releases.  Chemring  Ltd,  a  private  contractor,  operates  a  polar 
nephelometer  for  the  United  Kingdom  Ministry  of  Defense  that  has  been  used  to  measure  the 
phase  function  of  individual  millimeter  wave  chaff  particles  (1 5-/#m  diameter  by  4  to  8  mm  in 
length)  illuminated  by  a  94-GHz  collimated  beam.  This  system  is  capable  of  automatically 
acquiring  phase  function  data  for  a  wide  range  of  particle  orientations  relative  to  the  collimated 
incident  beam. 

The  Canadian  Defense  Research  Establishment,  Vaicartier  (DREV)  has  obtained  voluminous 
multispectral  extinction  data  in  a  large  laboratory  chamber,  and  extensive  near  infrared  lidar 
signal  returns  and  laser  transmittance  in  field  measurements.  Data  have  recently  been  reported 
that  correlated  concentration  measurements  with  integrated  back  scatter  returns. 

Data  acquired  by  the  Australian  Ministry  of  Defense  include  a  modest  number  of  chamber 
measurements  and  field  test  data  for  inventory  obscurants  extinction  and  for  graphite  fibers  used 
for  millimeter  wave  obscuration. 
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5.2  Scattering  Properties 


For  electromagnetic  wave  frequencies  of  the  order  of  100  GHz  and  lower  it  is  possible  to 
measure  single  particle  scattered  wave  amplitude  and  phase  as  a  function  of  particle  orientation 
relative  to  the  incident  beam  and  scatter  angle.  There  are  three  known  facilities  that  can 
accomplish  these  kinds  of  measurements: 

1.  The  microwave  scattering  facility  at  the  Institute  for  Space  Science  and  Technology, 

2.  The  bistatic  microwave  scattering  facility  at  the  Houston  Advanced  Research  Center,  and 

3.  The  United  Kingdom  95-GHz  polar  nephelometer  operated  by  Chemring  Ltd. 

The  ISST  microwave  scattering  facility  operates  at  a  wavelength  of  3.1835  cm  (9.4192  GHz) 
and  obtains  measurements  of  1)  extinction  (amplitude  and  phase)  (scattering  angle  is  0°),  2) 
angular  distribution  of  scattered  amplitude  (for  0°  elevation  and  0  to  172°  in  azimuth),  and 
3)  scattered  intensity  versus  target  azimuth  variation  for  fixed  scatter  angle  and  antenna 
polarization.  Figure  D-7  shows  the  ISST  microwave  scattering  facility  (Wang,  Microwave 
Scattering  Measurements  of  Small  Particles  of  Various  Target/Array  Shapes,  to  be  published). 
This  facility  was  designed  to  acquire  data  in  the  resonant  scattering  region  where  wavelength 
is  comparable  to  object  dimensions.  Thus,  data  typically  reported  from  this  facility  are  for  objects 
with  major  and  minor  dimensions  of  the  order  of  centimeters.  Original  applications  of  these  data 
were  tests  for  theoretical  models  of  irregular  particle  scattering  at  much  higher  frequencies. 
Scattering  object  indices  of  refraction  (predominantly  plastics)  were  measured  at  the  scattering 
facility  by  measuring  the  complex  dielectric  constant  using  the  standing  wave  technique  of 
Roberts  and  von  Hippel  (1946).  The  dielectric  constant  is  found  by  analyzing  observed  shifts 
in  standing-wave  maxima/minima  positions  as  well  as  the  voltage  standing  wave  ratio  patterns 
as  the  dielectric  sample  is  introduced  into  the  wave  guide-slotted  line.  Real  indices  of  refraction 
for  the  measurements  ranged  between  1.1  and  1.8.  Imaginary  components  of  the  scattering 
object  ranged  between  0.005  and  0.2.  Data  are  also  reported  for  aluminum  where  the  refractive 
index  is  reported  as  infinity. 

The  bistatic  radar  microwave  scattering  facility  at  the  Houston  Advanced  Research  Center 
(HARC)  operates  with  a  variety  of  polarized  microwave  sources  that  range  in  frequency  between 
2  and  40  GHz.  HARC  has  used  both  vector  network  analyzers  and  narrow  band  polarized 
beams.  Figure  D-8  shows  a  schematic  of  the  scan  geometry  of  the  HARC.  The  transmitter  and 
receiver  antennas  are  mounted  on  trolleys.  The  trolleys  can  move  along  their  respective  support 
trusses,  which  are  constant  radius  semicircles.  The  transmitter  truss  is  fixed  with  sufficient  offset 
to  allow  the  receiving  antennas  to  rotate  about  a  vertical  axis  between  0°  and  190°.  The 
transmitter  antennas  always  remain  at  0°  azimuth  throughout  their  entire  range  of  elevation.  The 
range  of  elevation  for  the  transmit/receive  trolleys  is  25°  to  150°.  The  scattering  object  pedestal 
mount  is  capable  of  rotating  the  object  through  400°  in  azimuth  with  an  uncertainty  of  ±  0.05°. 
HARC  was  designed  to  measure  bistatic  radar  cross-sections  of  objects  with  major/minor 
dimensions  of  the  order  of  centimeters  or  greater.  The  scattering  object  positioned  is  capable 
of  supporting  500  LB  radial  loads  and  2000  LB  axial  loads.  A  Ph.D,  dissertation  was  recently 
completed  on  a  12  component  error  model  for  HARC  (  Jersak,  1993).  Data  are  reported  for  a 
frequency  range  of  6  to14  GHz.  Jersak’s  work  represents  the  most  complete  known  error 
analyses  of  measurement  capability  for  any  of  the  facilities  described  in  this  report.  Figure  D-9 
illustrates  the  accuracy  and  precision  of  calibrated  object  measurements  at  HARC  and  compares 
the  results  with  theoretical  predictions.  Theory  and  experiment  agree  within  fractions  of  a 
percent. 
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Figure  D-7.  ISST  Microwave  Scattering  Facility. 


Figure  D-8.  Sketch  of  HARC  Bistatic  Microwave  Scattering  Facility. 
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Figure  D-9.  HARC  data  comparisons 


The  United  Kingdom  polar  nephelometer  operated  by  Chemring  Ltd  uses  the  highest 
electromagnetic  wave  frequency  (95  GHz)  of  the  three  facilities  described  in  this  section.  The 
polar  nephelometer  uses  Gaussian  beam  mode  optics  in  a  polarizing  Michelson  Interferometer 
to  obtain  plane  wave  illumination  of  the  scattering  object  placed  at  the  beam  waist.  The 
interferometer  geometry  allows  scattered  wave  phase  amplitude  in  the  fon/vard  direction  to  be 
measured.  These  data  are  used  to  compute  the  particle  extinction  cross-section  without  the  use 
of  calibration  targets.  The  largest  particle  that  can  be  examined  with  sufficiently  uniform 
illumination  is  approximately  22  mm.  The  lower  limit  in  size  is  defined  by  the  limiting  signal  to 
noise  ratio  which  is  determined  by  the  particle  support  system.  Past  measurements  have 
included  aluminum  coated  glass  fibers  approximately  4  mm  long  by  15//m  in  diameter.  The 
instrument  is  designed  to  automatically  scan  the  scattered  radiation  pattern  in  the  azimuth  plane 
and  to  position  the  particle  through  a  full  range  of  orientations  relative  to  the  incident  beam. 
Azimuth  scan  range  is  approximately  10°  to  170°.  Figure  D-10  shows  a  schematic  of  the  design 
of  the  instrument  (Poulson).  Figure  D-1 1  shows  examples  of  phase  functions  measured  with  the 
Chemring  instrument. 

All  the  nephelometers  described  in  this  section  are  capable  of  measuring  the  horizontal  and 
vertical  polarization  components  of  the  scattered  field.  It  is  only  necessary  to  measure  these  two 
components  to  fully  characterize  the  polarization  characteristics  of  the  scattered  field  when  both 
scattered  phase  and  amplitude  are  measured. 

5.3  Particle  Orientation  Distribution 

Nonspherical  particles  can  assume  a  preferred  orientation  after  dispersal  into  the  atmosphere. 
For  example,  it  is  theorized,  and  some  chaff  and  snow  data  show,  that  cylindrical  particles  orient 
in  the  atmosphere  so  that  the  long  axis  of  the  cylinder  is  parallel  to  the  ground.  Atmospheric 
turbulence  affects  the  particle  orientation  by  tilting  the  particle  long  axis  toward  the  vertical 
making  the  particles  fall  faster  and,  as  a  result,  attempt  to  reorient  parallel  to  the  ground.  This 
gives  rise  to  time  dependent  orientation  distributions.  For  orientation  distributions  to  occur,  the 
particles  do  not  have  to  be  cylindrically  shaped.  Even  raindrops  are  squeezed  into  oblate 
spheroid  shapes  due  to  wind  resistance  as  they  fall.  This  shape  gives  rise  to  droplet  orientation 
distributions  when  wind  cants  the  axis  of  symmetry. 

Scattering  from  irregularly  shaped  particles  depends  strongly  on  the  particle  orientation  relative 
to  the  polarization  vector  of  the  incident  beam.  Microwave  and  millimeter  wave  transmission 
through  preferentially  oriented,  naturally  occurring,  hydrometers  is  important  to  the  design  of 
communications  systems,  which  may  operate  in  two  characteristic  polarization  modes.  The 
polarization  modes  must  be  chosen  to  optimize  transmission  and  minimize  information  leakage 
between  the  modes.  The  polarization  mode  choice  strongly  depends  on  particle  shape  and 
orientation.  Similar  effects  and  problems  occur  for  the  design  of  millimeter  wave  and  microwave 
sensor  systems  operating  against  man-made  obscurants.  Therefore,  particle  orientation 
distributions  play  a  key  role  in  the  proper  design  and  application  of  millimeter  wave  and 
microwave  sensor  systems  in  the  atmosphere  and  from  satellites. 

Significant  research  programs  have  been  conducted  to  develop  theories  that  predict  particle 
orientation  distributions  for  raindrops.  However,  very  little  information  has  been  acquired  for 
orientation  distributions  of  obscurant  particles  in  the  atmosphere.  This  information  deficiency  is 
a  major  gap  in  data  needed  to  develop  and  verify  propagation  models.  Millimeter  wave  and 
microwave  system  designers  require  these  data  to  optimize  system  performance  against 
obscurant  countermeasures.  Obscurant  developers  need  these  data  to  develop  effective 
millimeter  wave  and  microwave  countermeasures. 
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Figure  D-10.  United  Kingdom  94  GHz  polar  nephelometer. 
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Figure  D-11.  Example  of  Chemring  polar  nephelometer  calibration  data. 


5.4  Particle  Number  Density  Size  Distributions  after  Atmospheric  Release 

To  estimate  transmittance  through  natural  hydrometers  or  manmade  obscurants,  it  is  necessary 
to  know  the  particle  number  density  size  distribution  and  the  numeric  spatial  distribution  along 
the  line  of  sight  between  transmitter  and  receiver.  Chamber  measurements  at  ERDEC  are  taken 
of  particle  dosages  (time  integrated  values  of  mass  concentration).  Dosage  is  used  to  infer 
concentration  in  the  chamber  by  assuming  a  homogeneous  distribution  in  the  chamber  during 
the  dosage  sample  time.  Chamber  measurements  are  also  made  of  concentration  as  a  function 
of  time  using  instruments  calibrated  to  respond  to  concentration  for  particular  obscurant  types. 
When  particle  size  and  density  are  known,  particle  number  density  can  be  estimated  from 
concentration  data.  When  concentration  is  sufficiently  low,  individual  particles  can  be  counted 
with  calibrated  light  scattering  instruments  to  measure  the  particle  number  density  directly  for 
monodisperse  millimeter  wave  obscurants  and  size  distributions  for  visual  and  infrared 
obscurants. 

5.5  Material  Dielectric  Constants 

The  complex  index  of  refraction  is  a  fundamental  input  required  by  virtually  all  electromagnetic 
wave  scattering  models.  Detailed  spectroscopic  data  determining  material  index  of  refraction  to 
a  wavelength  of  approximately  55  //m  are  available.  For  millimeter  wavelengths,  there  are  no 
direct  measurements  of  the  dielectric  constants  for  obscurant  materials.  In  this  case,  the 
complex  index  of  refraction  must  be  computed  from  the  conductivity.  High  frequency  dielectric 
constant  data  can  be  extrapolated  to  low  frequency  data  using  the  Drude  theory  of  metals.  In 
the  case  of  materials  with  a  permanent  dipole  moment,  a  modified  Debye  model  can  be  used 
for  extrapolation.  The  accuracy  of  these  extrapolations  is  not  known.  Data  are  available  of  fiber 
conductivity  at  millimeter  wavelengths.  However,  transformation  equations  must  be  used  to  infer 
complex  index  of  refraction  from  the  conductivity  measurements. 

Measurements  of  fiber  materiameter  wavelengths  are  of  materials  in  tow  (continuous)  form. 
These  measurements  are  made  in  a  resonant  cavity  wherein  the  cavity  Q  depends  on  the 
conductivity  of  the  fiber  material  pulled  through  it  (Goldberg).  Measurements  made  by  sampling 
long  lengths  of  tow  material  show  that  distributions  of  conductivity  exist  for  these  materials. 
Uncoated  materials  have  very  narrow  distributions  of  conductivity  values.  However,  coated 
materials  (for  example,  aluminum  coated  glass,  nickel  coated  graphite,  iron  coated  graphite)  yield 
broad  distributions  in  conductivity  values.  The  photographs  shown  in  figures  D-12  to  D-13 
illustrate  why  this  is  so.  The  figures  were  kindly  provided  by  Dr.  Ira  Goldberg,  Rockwell 
International  Corp.  The  figures  are  scanning  electron  microscope  photographs  obtained  by  Jane 
Hanamoto  of  Rockwell  of  obscurant  fiber  examples.  The  figures  show  conductive  layer  coated 
glass  and  graphite  fibers,  and  uncoated  graphite  fibers.  The  photographs  reveal  uneven  coating 
thicknesses  and  coating  breaks  and  cracks  along  the  length  of  the  fiber,  and  striations  in  the 
surface  of  the  uncoated  fibers. 

In  figure  D-12,  the  left  photograph  shows  cross-sections  of  typical  3//m  diameter  fibers.  The 
middle  photograph  shows  "dog-bone"  cross-sectional  shapes  from  Fortafil  Corp.  The  average 
diameter  is  6  to  7  //m.  The  right  side  photograph  shows  6  //m  diameter  graphite  fibers  coated 
with  nickel  produced  by  Americal  Cyanamid.  Note  the  variation  in  coating  thickness  which  is 
evident  even  on  this  scale. 
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Figure  D-12.  Three  microphotographs  of  graphite  fibers  (materials  discussed  at  the 
CIAOS II  workshop)  showing  end  on  views  of  fiber  candidates  for  millimeter  wave 
obscurants. 
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Figure  D-1 3.  Photomicrographs. 
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5.6  Obscurant  Characteristic  Poiarizations 

For  every  particle  ensemble,  there  exists  an  incident  wave  polarization  state  in  which  the 
polarization  is  unchanged  during  transmission  (Beckman,  1968).  The  incident  wave  polarization 
state  !3  called  the  characteristic  polarization.  Characteristic  polarization  is  a  useful  parameter 
for  defining  the  depolarization  of  electromagnetic  waves  transmitted  through  particular  types  of 
obscurants  and  the  atmosphere  in  which  they  are  dispersed. 

Characteristic  polarization  data  give  sensor  designers  guidelines  for  configuring  antennas  and 
polarization  geometry  to  optimize  sensor  performance,  the  obscurant  developer  information  on 
optimizing  obscurants  toward  one-way  obscuration,  and  obscurant  users  information  on  optimum 
atmospheric  conditions  for  effective  obscurant  performance.  Virtually  no  characteristic 
polarization  data  for  millimeter  wave  and  microwave  obscurants  exists.  It  is  also  important  to 
note  the  possibility  that  systems  optimized  for  obscurant  polarization  effects  may  be  inconsistent 
with  systems  optimized  for  natural  meteorological  effects. 

The  effects  of  particle  orientation  on  the  attenuation  of  polarized  transmission  in  the  millimeter 
wave  bands  is  limited.  Experimental  data  show  a  strong  polarization  dependent  effect  on  the 
mass  extinction  coefficient.  The  best  examples  of  these  effects  are  from  work  performed  in  the 
United  Kingdom  by  Edwards,  Poulson,  and  Butters.  Figure  D-14  plots  the  measured  mass 
extinction  coefficient  for  horizontally  and  vertically  polarized  35  GHz  radiation  attenuated  by 
aluminized  glass  cylinders  25  //m  diameter  by  3.1  mm  in  length.  These  data  were  meas  ured  in 
a  laboratory  facility  where  the  fibers  fell  from  rest  in  still  air  so  that  the  long  axis  of  the  fibers 
were  oriented  parallel  to  the  floor.  These  data,  plotted  as  a  function  of  mass  concentration  along 
the  transmissometer  line  of  sight,  also  show  potential  multiple  scatter  effects  because  of  the 
slight  decrease  in  mass  extinction  as  a  function  of  concentration.  Figures  D-15,  D-16,  and  D-17 
show  the  results  of  a  proof  of  concept  field  trial  for  25  /urn  diameter  aluminized  glass  fibers  of 
various  lengths  disseminated  from  a  RR170  cartridge  chaff  dispenser  for  aircraft.  The  1 .5-mm 
fibers  were  loosely  packed  providing  23  g  of  obscurant  per  cartridge.  The  4.0-mm  fibers  could 
be  loaded  as  aligned  packages  in  the  munition  to  provide  a  payload  of  76  g  per  cartridge  for 
aluminized  glass  and  59  g  for  nickel  coated  graphite.  The  transmittances  (shown  in  figures  D-1 5, 
D-1 6,  and  D-1 7)  were  produced  by  firing  six  cartridges  to  produce  an  obscurant  cloud  that  drifted 
through  a  transmitted  right  circularly  polarized  beam,  which  was  returned  to  the  source  using  a 
calibrated  trihedral  target.  Co-polar  and  cross-polar  returns  from  the  target  were  measured. 
These  data  show  typical  changes  in  signal  levels  for  a  two-way  transmission  through  the  cloud 
with  the  range  gate  of  the  sensor  locked  on  the  calibration  target.  The  aluminized  glass  and 
nickel  coated  graphite  fibers  affect  the  right  polarized  signal  power  differently.  The  aluminized 
glass  produced  an  increase  in  detected  signal  power  after  initially  decreasing  the  signal  power 
below  the  threshold.  However,  the  nickel  coated  graphite  only  attenuated  the  transmitted  energy. 
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Figure  D-14.  Laboratory  measured  35  GHz  mass  extinction  coefficient  as  a  function  of 
aerosol  mass  concentration  for  25-  ^m  diameter,  3.1-mm  long  aluminized  glass  fibers. 


102 


^  10 


$ 

^  -10 
C 

I 


Change  k* 


Time  (  seoonds  ) 

Figure  D-17.  Change  in  signal  return  from  a  calibrated  target  due  to  dissemination  of  4 
mm  long  coated  graphite  fibers. 
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AMSRL  BE  S 

BATTLEFIELD  ENVIR  DIR 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  BE  W 
BATTLEFIELD  ENVIR  DIR 
WSMR  NM  88002-5501 
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ARMY  RESEARCH  LABORATORY 
SURV  LETH  ANALYSIS  DIR 
AMSRL  SLE 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
SURV  LETH  ANALYSIS  DIR 
AMSRL  SL  S 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
SURV  LETH  ANALYSIS  DIR 
AMSRL  SL  C 
APGMD  21005-5066 

DEPT  OF  THE  AIR  FORCE 
7TH  WEATHER  SQUADRON 
APO  NY  09403 

USAF  ROME  LAB  TECH  LIB 
FL2810  CORRIDOR  W  SIT  262 
RLSUL  DOCUMENTS  LIBRARY 
26  ELEC  PKWY  BLDG  106 
GRIFFISS  AFB  NY  13441-4514 

ARMY  FIELD  ARTILLERY  SCHOOL 
ATSF  TSM  TA  SS 
ATTN  MR  TAYLOR 
FT  SILL  OK  73503-5600 

ARMY  DUGWAY  PROVING  GRND 
STEDP  MT  DA  M 
ATTN  MR  CARLSON 
DUGWAY  UT  84022 

ARMY  DUGWAY  PROVING  GRND 
STEDP  MT  DA  L 
DUGWAY  UT  84022-5000 

ARMY  DUGWAY  PROVING  GRND 
STEDP  MT  M  L 
ATTN  MR  BOWERS 
DUGWAY  UT  84022-5000 


no 


DEFENSE  TECH  INFO  CTR 
8725  JOHN  J  KINGMAN  RD 
SUITE  0944 

FT  BELVOIR  VA  22060-6218 

OPTEC  TECHNICAL  LIBRARY 
4501  FORD  AVENUE  SUITE  820 
ALEXANDRIA  VA  22302-1458 

ARMY  FRGN  SCI  &  TECH  CTR 
CM 

220  7TH  STREET  NE 
CHARLOTTESVILLE  VA  22901-5396 

NAVAL  SURFACE  WEAPONS  CTR 
CODE  G65 

DAHLGREN  VA  22448-5000 

NAVAL  SURFACE  WEAPONS  CTR 
CODE  J41 

ATTN  MR  C  RUSSELL 
DAHLGREN  VA  22448-5000 

LOGISTICS  CENTER 
ATCL  CE 

FT  LEE  VA  23801-6000 
TACDOWP 

LANGLEY  AFB  VA  23665-5524 

ARMY  NUCLEAR  &  CHEM  AGENCY 
MONA  ZB  BLDG  2073 
SPRINGFIELD  VA  22150-3198 


Record  Copy 


TOTAL 


